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CHAPTER 1 


INTRODHCTIOH 

The East Eourier Traiisforn is a netliod for efficiently 
conputing tlie discrete fourier transforri of a sequence of 
data samples. This technique greatly reduces the numher of 
computations required to calculate such s. transform on a 
digital computer. Consequently;, it has made feasible the 
use of Eourier transforms in the analysis of many problems . 
tha,t were previously approached by other methods. Eourier 
transforms are now routinely used in sucii diverse areas as 
ceisij'iiic exploration, speech analysis, eclio^ranging systems , 
vibration analysis, image processing, and many others. 

Recent years have witnessed a great increase in the 
availabilitjr of small, relatively inexpensive computer 
(mini/micro-computer ) , which can be employed for East Eourier 
Transform of signals in real-time applications. Such 
computers may be used to sample incoming signals and to 
perform certain calculations using these data so that the 
results become known as the process continues. The speed 
of these computers may limit the number of sample points on 
which the PET can be done in a specified time. This restri- 
ction can be overcome either by interfacing the special 
hardware with these computers or by the use of special 
techniques like parallel processing, memory organisation, etc. 
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Thus zeeping in view the capability and in expensiveness 
of nicro-cornputers, it was projected to develop software 
and Hardware for 1024 points PPT using b . microprocessor. 

The i'TT package was designed specially to analyse data 
from a digital correlator for the Bake Troposca^tter System, 

At the output of digital correlator, 10 samples (10 Hormal 
and 10 Quadrature) are obtained in 6.2 msecs and as such 
1000 samples are accumulated per file in 620 msecs. Thus 
the time constraint of 620 msecs for IPT of a ^ile was 
fixed to achieve real-time environment. Out of this total 
time, most of the time is spent, in coraputo.tion (very small 
time is spent in house keeping, 1/ 0 operations and recording 
of the Fourier coefficients), which consists of addition and 
multiplication time. Since the time required for multiplica- 
tion is assumed to be much greater than tha.t for addition, 
the total time for computation is approximately proportional 
to the multiplico,tion time for the system. 

The available microprocessor either do not have hardware 
multiply instruction (Intel 8080, Motorala 6800), or the 
one like TMS 9900 which have it, perform the multiplication 
slowly. Thus it is essential to have a hardware multiplier 
unit interfaced with the microprocessor to achieve the 
time constraint of real-time applications. 
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There are four different ways ; (a) Sequential 
processor, (b) Cascade processor, (c) Parallel iterative 
processor and (d) Array analyzer, in which PFT processor can 
he organized [ 5 ]. The sequential processor is characteri~ 
zed hy one arithmetic unit and a computation tine proportional 
to ^ log2^j where h is the total number of samples. The 
cascade processor has log2l'T arithmetic units, and the compu- 
tation time is proportional to ^ in this case. The parallel 
iterative processor is characterized by arithmetic units and 
a computation time proportional to Log2l'^* Array analyzer 

is considered in which all ^ Log2H' operations are performed 
in parallel and the execution time is simply the time 
required for performing one basic operation. It ha,s 

TT 

2* log2h’ arithmetic units. 

Evidently out of all these different organization 
schemes, the sequential processor is the slow'est and at the 
same tine the cheapest one, while the Array anaJ.yzer is the 
fastest as well as the costliest one. The execu.tion time 
has been reduced in the other types of non sequ.ential proce- 
ssors by introducing parallelism in arithmetic operation. 

Since the sequential processor scheme is simple in its 
organisation and requires only one microprocessor, we 
considered to start with this organization. To counteract 
the inherent slowness of this scheme and thus to achieve the 
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time constraint for real-time application^ a compler liard-ware 
multiplier unit (having four multipliers in- it ) instead of 
simple multiplier unit is required to he interfaced with 
the microprocessor. 

Out of the two microprocessors (Intel 8080- and Motorala 
6800) available,, Intel 8080 was selected because of its 
better suitability for scientific calculations and available 
software support like cross-assembler and simulator on 
7044 computer. 

In this thesis, procedures are developed for implemen- 
ting Fast Fourier Transform on 8080 microprocessor. Special 
algorithms are devised that cut down the time required to 
calculate FFT, thus making them useful for real-time 
applications. 

A companion thesis by Mr. U.S. Bhaliat discusses the 
hardware related to the implementation of FFT on microproce- 
ssor [7]. 

The details of FFT algorithm and their basic properties 
are described in the second chapter. Third chapter contains 
the design of FFT program and its optimization. FFT imple- 
mentation on microprocessor and microprocessor limitations 
for real-time application are discussed in the fourth 
chapter. Conclusions and suggested future work are con- 
tained in the fifth chapter. 
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OHAPTER 2 


EPT ALGORITHI© 

PET is an algoritiim (i.o, a particvilar method of 
performing a series of computations) tlie^t can compute the 
ais Crete Pourier transform much more rapidly than other 
available algorithms. Brute force calculations required 
opera,tions since R Pourier coefficients are to be evsl-uated 
and each is the sum of If products. 

An algorithm developed by J.¥. Cooley and J.¥. Tukey 
reduces the computational load to IT log-glT ■where B is the 
base (typically a power of 2 such as 2,4,8 or 16) to which 
the logarithm of IT is taken and also represents the n'umber 
of data from the full set of hi which are processed in each 
substep of the procedure. 

There have been some variations (like S ond e-Tulrey 
algorithm) to the original algoritlim and at present there 
are many algorithms available for calculating the PPT, ¥e 
will inspect some of the variations of the basic PPT algori- 
thm and computational structxnro. It is stressed that 
algorithms presented is not intended to be a complete set of 
such algorithms. Indeed, there have been many additional 
modification depending on the particular requirements, the 
limitations of available hardware, and the ingenuity of 
individual designer or programmer. 
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Most of those algorithms may bo olassifiod as either 
(a) in-place or (b) natural input -output [2], In this 
chapter, these algorithms will be studied by the use of 
signal flow graph and their basic properties will be iden- 
tified, 

2.1 In-place Algorithm 

An in-place algorithm is one in which a given component 
of any intermediate vector may be stored in the same location 
occupied by the corresponding component of the preceding 
vector. This type of algorithm requires less total storage, 
but at the same time the computational time is higher than 
that required for natural input -output algox’ithm. Other 
peculiar characteristic of these algorithms is either the 
output spectrum appears in an ujina.tural order or they 
require that the input data be arranged before entering the 
computation array. Thus in-place algorithms require the 
reordering of input or output data. This reordering process 
is referred as scrambling operation and is discuss od next, 

2.2 Scrambling Operation 

The scrambled value of a given integer m wi.ll be defined 
as ms [2], Assume that m can be represented in binary form 
as 

“ ~ “^N-2 °^l^o 
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The scrni-biod value of m is defined as 

nig = mQ nil ^if-i 

Thus, file scrani’bled. value of a given infeger is 
a new numher obtained by reversing the order of all the bits 
in the binary representation of the given number. Note that 
if m is scrambled twice, the origina.1 value is obtained. 

Por illustra.tion, values of m and ms for If = 8 and 
If = 16 are given in decimal and binary form in Table 2,1. 

With some of the in-place algorithms, the data must 
either be scrambled before or after processing. If we 
assume that input is in natural order, then the output 
requires scrambling. Then, at a particular location m, the 
component appearing at the output is not }[(m), but rather 
X(ms). In this case, it would be necessaiy to go to location 
ms to obtain the component desired for the index m, 

2,3 Signal Plow G-raph 

PPT signal flow graph for in-place algorithm for If = 8 
is shown in Pig. 2.1. It consists of data array and compu- 
tational arrays. Data vector or array is represented by a 
vertical column of nodes on the left of the graph. The 
vertical columns to the right of data array correspond to 
computational arrays and in general there will be r 
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computational arrays whore r = 10^23’. 

2.4 Basic Properties of In-place Algorithm. 

Prom the flow graph, the following basic properties 
[4] of in-plo.CG algorithm can bo identified ; 

i) Number of computational arra.ys 

r = log 2 N 

ii) In each computational array every node has tvro in- 
coming paths and two outgoing paths 

iii) In each array, there are two nodes whose input pa,ths 
originate from the sojne pair of nodes in the previous 
array. Two such nodes are grouped as ’a dual node 
p8,ir' . In 8.ny array there are N/2 pairs of dual nodes. 

iv) Each array requires N/2 comples multiplications and N 
complex additions. Hence the total number of multiplica- 
tions and additions required are (N/2) log 2 N and 

N log 2 N respectively. The ratio of direct to PET 
computation time is 

N^ ^ 2N 

1/2 Log2N iog2N 

v) The computation of a dual node pair requires only one ■ 

multiplication and two additions. If the weighting factor 
at one of the nodes in a dual node pair is then the 

weighting factor at the other node of the pair is 

fP 



C,o'h7pj{ji-ta.i:t07iai 

J=z 1=3 



B act! Ti ode^ ”. X^CK) a^vcl + 

Ff T- — ^7^-P^a^e 
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# + lf/2 ^ 

Then 

Xi(lc) = + Er/2^) 

2:^(k + E/2^) = x^„i(k) - x^__^ (k + IT/2^) 

iiere x^(k) indicates k coniponent in tlie l”*"^ 8,rray. 

vi) The spacing hotween dual node pairs differ from ax’ras?' 

th -■ 

to array. In the £, array (v, = 1,2, r) the 

f P 

spacing is N/2' , i*e. x (k) and x (k + IT/2'' ), 

cjonstitute a dual node pa,ir. 

vii) To evaluate the value of P, the exponent of ¥, for any 
index in a given array, the following procedure is 
followed. Represent k, the node index in the >£th 
array, in binary forra td.th r bits, retain the most 
significant 'C bits and add (r ) loading zeros 

to form a r bit binary number. PLevorse the bit 
order of the resulting number. The decimal equiva3.ent 
of the final binary nnuber gives the index P. The 
weighting factor for the kth node of the Ith array 
is ¥^. ■ 

viii) The output after r arrays is in scrambled form. To 

unscramble the output x(k), x^rite the index k in binary 
form with r bits and reverse the bit order. The 
resulting decimal number is the index n of x(n). 



12 


2.5 Natural Input-Outiout Al-rorithn 

Natural iiiput-output algor5,tlim is one in whida a 
given component of any intermediate vector ma,y not be 
stored in the same location occupiod bj?" the corresponding 
component of the preceeding vector , thus requlririfg the 
extra memory for storing the intermediate result. An N 
point FFT (N real and N iraagiriary ) wilD. roqulre 4N uords 
of DiGiaory as compared to 2N words requ.ired for in-place 
algorithm. These algorithms, of course, nmintain the natural 
input-output order and thus do not require scrambling/ 
unscrambling at the input /output levels, and as such a.re 
faster as compared to the in-placo algorithm. 

2.6 Signal flow Graph 

,i signal flow graph for N = 8 for natural input-output 
algor ithra is given in Pig. 2.2. There are four columns and 
each column contain eight entries (Number of sample points). 
The variable (k) is used to denote the value of a given 
node in the array, where L is the number of column and k is 
the number of the component within the column. In genoral 
■I varies over the range 0 1- t 1 r with = 0 at the left 

and r = loggN, and k varies over range 0 l.k I N-1 with 
k = 0 at the top. 
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2.7 Basic Properties of Fatr.ral Inpu.t“Outpiit Algorithm 

Properties (i) through (iv) are same as for in-place 
algorithms [4]. 

v) The computation of a d^^al node pair requires oiil37' one 
multiplication and two additions. If the weighting 
factor at one of the nodes in a dual node pair is 

then the weighting factor at the other node of the pa.ir 

is wP + 

WP + Pf/2 = -V(P 

Then 

P 

x^'(k) = ^ IT/2”’ ) 

I 

X/ (k + lT/2) = (i + lT/2 ) 

i depends upon location of k in the arra 3 r and the numhor 
of array . ' 

vi) The spa,cing between dual node pairs is same for all 
arrajT’S and is 11/2. 

vii) To evaluate the value of P, the exponent of ¥, for any 
index in a given array, the following procedure is 
followed. Represent K, the node index in the- /th array, 
in binary form with r bits,- retain the most significant 

/ bits and add (r -1 ) zeros in front of it to form 
a r bit bimry number. The decimal oquivalont of the 
binary number gives the index P. The we?-ghcing fo.ccor 
for the kth node of the .-^.th array is ¥^. 
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viii) The output after r arrows is in natural order and 
thus does not require uns or unbling . 
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CHAPTIR 3 
P.FT PROGRA.M DESI® 

Programming reciuires a disciplined a,pproacii to the 
translation of requirements into unai'nbiguous instructions 
for a suitable computer. However, programming involves much 
more than merely transcribing some symbols, it co}:]Sists of a-t. 
least five major steps ; 

i) Design 

ii ) Coding 

iii) Translation 

iv) Testing 

v) Debugging 

In general, if the design is not done carefully, testing 
and debugging will take an inordinate long time to complete. 

The design of a computer program requires not only the under” 
standing of the problem, but also the suitable selection of 
algorithms. The selection of algorithm becomes quite critical 
in real-time environment, where not only the program should 
work , but should work efficiently i.e. should talce minimum 
time for its execution. This requires the optimization of the 
program before it is tested and debugged. The optimization 
puts a great strain on the programmer, because he has to 
consider the efficiency and under standibility of the program 


at the same time. 
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This chapter will discuss the design, of PET program 
and its optimization for its application in real-time. 

3.1 Analysis of Available EFT Program 

Before starting the design of our own EFT program, it 
was decided to analyse the available EFT programs. After 
surveying the literature,, only one program written in EORTRAM’ 
was found in the book 'The East Fourier Transform' by 
Brigham [1]. This program is given in Appendix 'A'. The 
program was run on 7044 compute for different number of sample 
points (1024, 2048, 4096) and execution times corresponding 

to these samples were calculated, by the use of Function Time 
(IDUM) available in FORTRAN. The results are given in 
Table 3.1. 

It can be seen from the ta.ble that the program takes 
24000 msec ', for 1024 points FFT. It is not worth trying 
to code this program into 8080 assembly language, since it 
will not be possible to achieve the time constraint of ^ 

620 msecs, for 1024 samples, nowever hard one may try to 
optimize this program. 

The slowness of this program is attributed to the 
following reasons i 

(a) The program uses in place algorithm which is slower 
than natural input-output algorithm. 
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(Id) -The weights ¥ (where ¥ = are oalculated by the 

use of library function £I¥S and COSH®, which consumes 
time and makes the program slow. 


ho. of Sample 

Execution Time 

N 

in msecs 

1024 

24000 

4096 

120000 

8192 

225000 


Execution^j^Time ~ EORTEAN Program 
(Ingham's Book) 

Table 3.1 

3.2 Selection of Algorithm 

As stated above, the in-place algorithm is slow and is 
thus not fit for real-time environment. As sudh, the natural 
input-output algorithm is selected for the development of 
FPT pro gram , 

3.3 EORTRAH Program Design 

The flow-chart for PET program utilizing natural input- 
output algorithm is given in Fig. 3.1. FORTRAN program is 
developed based on this flow-chart a,nd is given in Appendix 


? o ! 
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Input 

XREAL(I) 

:XBIAG(N) 


I1=LV2, IB=1 


I = 1,L 




E=1,IB 


10=0, ID=IA 

4 


M=ia,iD-i 


C=COS(lC),S=SIiT(IC) 




AREAL=XifflAL ( M+ 10) 

A IMA &=Xim 0 ( M+ 10) 

BREAL=XRBA1(M+ID) * 0 + Xim0(M4-ID)*S 
BB'IAMIMAtKMflD) * 0 - ZEBAI(M+ID)*S 
YEEAI(M) = AEEAL + BBEAl 
YIMAO(M) = AIMAG + BII-IAG 
YEaAL(M+M/2) = AHEAB - BjIEAL 
YBUG(M/2) = AIMAO - BIMAG 


IO=ID, rD=ID+IA 

4^ 


IA=IA/2, IB=IB3£2 
M=0, H-l 




XEEAL(M) = YREAL(M) 
XIMAG(M) = YIMAG(M) 






OUTPUT 

xbbab(u) 

XIMAUCl) 


Pig. 3,1 Plow diart - PPT natural input -output 
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As -^he program is based on natural input-output 
algorithm, it requires more memory for storing intermediate 
rt;Sults, lor 1024 samples (1024 rea,! and 1024 ima,ginary)j 
it requires 4096 words of data sticrrage as against 2048 words 
in case of in-place algorithm. YEEAl and YIMAO are used for 
storing the intermediate results. After each iteration, 

YREAl and YIMO are transferred to ZREAL and ZIMAO and 
process continues for r iterations (r = Iog 2 l)- finally, 
the results are stored in lEEAI and XIMG and the data 
which was initially in XREAl and ZII-IAG is lost. 

The program was run on 7044 computer for 1024 sample 
points and it talces 3300 msec, for its execution. This time 
is 14 percent of the time taken by Brigham’s program. The 
selection of natural input-output oYgorithm ha.s reduced the 
time from 24000 msecs to 3300 msec for 1024 points FPT 
(86 percent reduction). 

3.4 Weights of ¥ 

It is not required to call library function SHE and 
COSIHE to calculate the weights of ¥ (as has been done in 
Brigham's program), in each iteration. Yalues of SIKE and 
OOSIEE can be generated and stored in a array before entering 
the EFT. E/2 values of SI IE and H/2 values of COSINE are 
required to be stored from 0 to ^ (N/2 - 1) in a step 
size of 2it/E. This will correspond to E/2 weights of ¥, i.e. 
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In the program, instead of 

calling the library function, the index for SIME and aoSIIIB 
array is calculated and appropriate values of SIBTE and 
COSIITE are fetched from the array for Cralcula/bion of FET. 
This requires extra 11 words for storing the values of SIEE 
and GOSIEE, 

Latter, in the course of development of p^eogram it was 
found that ¥/2 values of SIEE and IT/2 values of COSINE are 
not required to be stored. Only lT/4 1 values of SINE from 
0 to Tc/ 2 in a step size of Oii/ll are needed and the other 
values of SINE and all values of C0SIII3 can bo generated from 
these values. These N/4 + 1 values of SIITS correspond to 
N/4 + 1 weights of ¥, i.c. ¥°, ..... 

For example, for N = 1024, 257 values of SINS are stored 
from 0 tc %/2 in step size of 2'n;/l024 shown in Fig. 5.2. 

rm T _L T tO t t^L T T ^ 5 S 

These coi're spend to ¥ , ¥ 4 

¥e will elaborate how values of SINE and COSINE are 
ca,lculated from the values of SINE stored. 

5.4.1 Calculation of SINE values (N = 1024) 

For calculating the values of SINE required in the FPT 
program, following procedure is adopted : 

i) The program should check whether the index is 1= 256 or 

> 256 . 



Location in 
memory 


Ya.lue of nine 
stored 


0 0 

1 2t:I102A ± 1 

2 2n;/l024 3: 2 


256 


2%/l02A * 256 

(n/2) 


Values of SINE Stored in Memory 
Figure 3.2 



ii) If index frr 256, then rise the index aa it is to 

fetch the appropriate valu-c of from tho array. 

iii) If index 256, the index is subtracted from 514 

(in general from 2(11/4 + 1)) said this vaJ.ue is used as 
the index for fetching the appropriate vrJ.ue of SINE, 

3.4.2 Calculation of COSINE Values (N = 1024) 

The following procedure is followed ; 

i) The program checks whether the value of index 256 
or > 256 (in general, index N/4 or .> N/4). 

ii) If index 256, subtract the index from 256 (in 
general, from 11/4) and use this value as index to fetch 
the value from SINE array. 

iii) If index > 256, subtract 256 (in general N/4) from 
the index and use this ve,lue as tho index to fetch 
the value from SINE arrry. This will give appropriate 
COSINE values. 

This method of storing only (N/4 1) values of SINE 

instead of storing N/2 values of SINE and N/2 values of COSINE 
will cut down the memory requirement for storage of weights 
of ¥ from N words to N/4 + 1 words. 

The program given in Appendix 'B' was modified to 
incorporate the idea of weights storage xnstead oi calling 
library function to calculate them. This program is given in 
Apnendix ’ C , The program vras run for 1024 sojaplss and the 
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oxGCution '’jime was oalculated. The program takes 2450 msoo 
and this technique resulted in 26 percent reduction in 
exocLition time. 


3.5 Specia.1 loops 

The program developed was further optimized by exploiting 
the very nature of the algorithm. There are r iterations 
(r -- log 2 ir) in the FTT program. In the first iteration, 
the wcigiit of \'J is zero, which gives value of Sllffil as zero 
and GOSIFE as one. This means that the components in first 
array can be computed by simply addition and subtraction and 
the no id for multiplication with SIFE and GOSIKB is eliminated, 
For cx.ample i 

%^^^^imag “ ^o^^'^imag ^o^"~ ^-^^^imag 

X^Ck + F/2)^^^^ = ^o^l^^roal " ^o^^^ ^ 


X^(k + ^o^^^^imag 


Xo(k 


Similarly in second iteration, the weights of ¥ are 0 and 
F/4. Zero weight of ¥ gives values of S'lHE and GOSIFE aS 0 
and 1 respectively, whereas F/4 power of ¥ will . give value of 
SIFE as 1 and GOSIFE as 0. As such, the components of second 
a,rray can also bo calculated without any multiplication. 
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The program was then modifiod such that components in 
array 1 (iteration 1) and array 2 (iteration 2) can be 
computed separately before entering the main FPT programs, 
which will now be executed for (r - 2) number oftimes. 

The modified program is given in Appendix 'D' , This 
technique reduced execution time for 1024 points PfT from 
2450 msecs to 1900 msecs. 

3.6 Memory Swapping 

In natural input-output, the components calculated 
during the iteration are stored in the intermediate locations. 
Those components are transferred back to their original 
locations before starting the next iteration. In the program 
Y.RBAL and YIMAG are used as intermediate locations and 
ZIlEAl and ZIMAG are the original locations. As such 
2N (F real and W imag) memory transfers are required after 
each iteration and thus a total. of IHr (r = Log 2 N) memory 
transfers for PPT program. 

These memory transfers can be avoided by the memory 
swapping technique. In this technique, one works on 
(XREAL, XII-IAG) and stores the result in (YEEAL, YIMAG) in 
the first iteration. How instead of transferring YRnAL and 
YIEAG to XREAL and XIIIAG respectively before starting the 
next iteration, one has to work on (YRBAL, YIMAG) and store 
the result in (XREAI, XIMAG) in 2nd iteration. In the next 
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iteration, one has to work on {iKEAJj, XIII/iG) ojid store the 
result in (XRBAL, YIII/VG-) and so on. Fina.ll}/', the result will 
be in (TIffiAL, XBUG) if r is even .and in (YEEAl, YimO-) if 
r is odd. 

This technique, in case of r being even, will out down 
the memory transfers from 2Kr to zero and in case of r being 
odd, to 2N transfers only. 

Tho technique may not save much time for lar-ge computers, 
in which memory transfers are quite fa,st, but for mini/micro 
computers and especially for micros, which do not have 
memory to memory transfer instruction, the saving in time 
will be appreciable. 

The progra.m was then modified to utilize this technique 
and resulted in execution time reduction from 1900 msecs 
to 1500 ms CCS . 

The benefit gained from tho use of techniques discussed 
3.2 tlirough 3.6 are shown in Table 3.2. 

Thus the selection of algorithm, exploiting the very 
nature of algorithm emd some ingenuity helped in the 
development of optimized program. 
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Tocliniquo 

Basic Di’isham Progrojj 

Selection of algoritlini 
(from in-place to 
no.tirral ) 

Weight storage 

SpociaJ. loops 

itinory swapping 


Exocution Tim.e Time- 
for 1024 points 

24000 rasecs 100 

3300 msecs 14 

2450 ElSOCB 10 

1900 msecs 8 

1500 msecs 6 


Execution Time Reduction 


Table 5.2 


reduced to 
__ ^ 

(rofereiice) 
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aapTER 4 

MICROPROCESSOR IMPIEMTATIOH 

Ev. r GincG the development of the first microprocessor 
in 1971, there has been a trenendous increase in the appli- 
cation of microprocessor in diverse areas such as process 
control, instrumentation, consumer products, da,ta acquisition 
and many rcal-tino application. In most practical cases, 
it is realistic to consider /ap for any hard wired logic 
cmployiiig more than 50 or 60 ICs, having more than a trivial 
number of steps in the flow chart and having some logical 
and aritiimetic data processing requirement, 

Eventhough microprocessors have got an advantage in 
their usage as compared to hard wired logic, they have their 
limitations in scientific and roal-tirae application. These 
limitations mainly arise because of the word length and the 
speed of available microprocessors, which in turn dictate 
the available accuracy and the very utility of microprocessors 
in real-time application. 

In this chapter, we will discuss implei-ientation of EET 
on 8080 microprocessor and microprocessor limitations for 
scientific real-time application, because of word length 
and speed. 
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4.1 Sclv.’ction of Pixod-Point Arithmotic 

To iinplonent the PPT on. 8080 rnicroprocossor , wo sele- 
ctod fixed-point nrithmotic over floating point ai'itlimetic, 
because fixed point is faster than floating point though less 
accura-to. In tlio PPT program utilizing fixed point exithme- 
tic, tlio input sequence (data) is scaled such that it can be 
reproBonted by B bits plus sign and the binary point is assumed 
to lie to the loft of the loftraost .magnitude bit. As we 
novo from stage to stage of the program, the magnitudes of 
the numbers in the sequence goner o.lly increase which means 
that there is possibility of incurring overflows dui-ing 
difforont stages of computations. To prevent overflows, some 
technique of scaling is required in fixed point arithmetic. 

4.2 Scaling 

The inner loop of the power of two P.PT algorithm 
operates on two complex numbers from the seq^ionce [6]. It 
tiikos these two nuribers mad produces two new complex nmabers 
which replace the original ones in the sequveiaoe, let 2^(i) 
and X^(o) "be the original complex numbers. Then, the new 
pair are given by 

\+lU) = 
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At r-ich stagG, tlio algorithm goes through the entire 
saqu.nico of N niurfoors in this fashion, two at a time. If 
n = 2 % tiion tho nuraher of such stages in the coriputation is 


With tho assumption that the binary point lies at the 
ostreno loft, the re.lationship om ong tho numbors in mth 
stage end n+lst stage- is as shown in lig. 4.1. fhc outside 
square gives tho rogioii of possible values, Rc [S^(i)] <1 
and circle inscribed in this square gives 

the rogionj X^(i) | 1. Iho inside square gives tho region 
Rc[Xj^(i)] < V* finally, the circle inscribed 
in this latter square gives tho region 

(i) and inside the smaller circle, then 

Xj^_i, 2 (i) and will bo inside the larger circle and 

hence not result in an overflow. Consequently, if wo 
control the sequence at the mth stage so that \j(i) 
wo are cortain wo will have no overflow at the m-i-lst stage, 
ITowevor, if Xj.^(i) and X|^(o) are inside the smaller square, 
then it is possible for \+i(i) ) outside 

tho large square and hence result in overflow. Oor^equently , 
WG can not control tho sequenco to prevent overflow by 
keeping the absolute values of the- real and imaginary parts 
less than one-half. 





Re ICL-fc LO^'iShip Bei Nombe.^/ Iv 7>UK 
CLVcl 1'T\b±s-t e 


ptqtL)^e 4-1 


1 I.T. KANPUR 
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Tho three techniques of scaling, which when applied 
prevent overflow, are given below : 

i) Shifting Right one Bit at Every Iteration 


If the initial sequence 1^(1 ), is scaled so that 
i X (i) !< 4 for all i and if there is a right shift of 
one bit after every iteration, then there will be no 
overflow. i 


ii) Controlling tho Sequence so that 


\i±) k i 


Aga,in assurae tho initial sequence is scel.ed so that 


"0 
clicdc 


X^(i) I < -!■ for all i. Then at each iteration we 

Xj.^(i)l and if it is greater than one-half for 
any i we shift right ono bit before calculation through- 
out the next iteration. 


iii) Testing for an Overflow 

In this case the initial sequence is sccRed so that 


RG[XQ(i)] < 1 and IJX^(i)] < 1. Wlienevor an overflow 
occurs in an iteration the entire sequence (part of 
which will be new results, part of which will be entries 


yet to be processed) is shifted right by one bit and 
the iteration is continued at the point at which the 
overflow occured. 


The first technique is the simplest and easy to adapt 
for microprocessor. This method gives less accuracy than the 
other two techniques, since it is not generally necessary to 
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rescale the sequence at each iteration, there is an un- 
iiccessaDry loss in accuracy. The second technique requires 
chocking during eo.ch iteration and will take good 

aiiount of tine if inploraented on microprocessor. Since we 
had the tine constraint, we decided not to go for this 
nothod. The third technique requires checking the overflow 
and as such the machine selected should havo overflow as 
ono of its status flags for its implenentation. Since the 
8080 microprocessor docs not have an overflow flag, this 
toclniiquo also could not be adoptod. As such, we selected 
the first technique of scaling to bo u.tilized in the FPT 
program to be implemented on 8030 microprocessor (o.ccuracy 
sacrificed for the sake of saving in execution tine). 

4.5 8080 Assembly Isnguago Coding of PFT Program 

The program designed in Chap tor 5 was coded into 
assembly language of 8080 microprocessor. Since the data 
from digital correlator is 5 bits plus sign, wo assumed that 
8-bits word length of 8080 micropx'o cossor mil be enough. 

At this stage, we hid not analyse whether 8 bits (7 bits for 
magnitude and ono bit for sign) word length will be enough 
from accuracy point of view, later simulation proved that 
8-bits word length is not enough. 

4.3.1 Precautions for assembly language coding 

The program has got three loops a) outer loop, 
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b) inner loop and c) inner 17080 loop as slioirn in Fig. 4.2. 

Tlie instructions in the outer loop get executed by 
r (r = Log2^0 ntiraber of tines. Tliat is if IT = 1024, the 
instructions in this loop will be executed by 10 tiraes. 
is such this loop is not at all critical fron selection of 
instructions to be used in this loop. 

The inner loop gets executed for (IT— 1) tines and requires 
none precaution in the selection of instructiouis which take 
less tine for their execution. 

Since there anre no coraputation involved in outer and 
inner loops, the tine taken for their execution will be a 
small proportion of the total execution tine of FFT program. 

The most critical of the three loops is the inner most 
loop. The coraputation of the members of the array, two at 
a time for power of two algorithm, is carried out in this 
loop. For the calculations of the entire members of the 
array, this loop is executed lT/2 tim.es. Since there are r 
arrays in the FFT program, this loop will be executed 
(F/2 r) times. For IT = 1024, this loop will be executed 
5120 number of times. The execution tine of this loop will 
be more or less equal to the execution time of the FFT pro- 
gram. Every care should be taken in the selection of 
instructions to be used in this loop. 
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4.3.2 Ro-u^h Bstiiiiatioii of Execution Time 

The execution time of eacli instruction (time taken 
from 8080 UrJER' s MA.ETJAL) was added in the outer, inner and 
inner most loops. The time calculated for each loop was 
multiplied by the number of times that loop gets executed 
and the total time, which is the sum of times taken by 
tJiree loops, gave us the rough estimation of the execution 
time for PET program. The tine estimated was 6000 msecs. 

This value is a lower bound, because house keeping and 
input“Output operations have not been included in this 
estimate. 

4.3.3 Interfacing Hardware Multiplier 

As stated above, the rough estimation of execution 
time gave us 6000 msecs. This time is vjith software multiply 
(8 X 8) routine, which talces 250 psecs for 8x8 bits 
multiplication. It is clear that tine target of 620 msecs for 
1024 points PET can not be achieved with software multip3.y 
routine. As such, it is required that EFT prograia should 
have hardware multiplier unit to meet the time constraint. 

4.3.4 Time Estimation with Hardware Multiply 

The estimation of execution time of EFT program with 
Hardware multiply was made. Here we assumed that a complex 
multiplier unit, which have 4 multiplier units, is inter- 
faced. This multiplier can do 4 multiplications at the same 
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■fciue and .'■'iveG back two rosul-ts after multipbi cation. ¥e 
cad. culaced that a miiiini'ujTi of 30 )..isecs will be needed to 
perforu the above operation, because 4 'OUT’ instructions 
will be needed to send 4 operands aud 2 ' lU' instruction for 
receivers the results (each OUT/'II instruction takes 
5 jiseco). The execution time with the complex hardware multi- 
plier was calculated as 1000 msecs. 

The 1000 msecs time calculation is with the available 
8080 A microprocessor which has a clock frequency of 2.08 MHz. 
Other faster version of 8080 microprocessor, which have clock 
of 3.00 MHz, will give around 70 percent of the estimated 
time. As such we estimated that the time target of 620 msecs 
required for real-time application of fFT can be achieved with 
faster version of 8080 microprocessor interfaced with complex 
hardware multiplier unit and vrith a little optimization of the 
developed program. 

4.4 Simulation of 8-bit Machine (microprocessor) on IBM 7044 

Before going ahead with the testing of assembly language 
program, written using 8-bit data length, we decided to 
check whether 8-bits word length will be enough from accuracy 
point of view. ■ 

A I'OHTIIAI program which simulates the 8-bit microprocessor 
on IBM 7044 computer and incorporates fixed-point arithmetic 
was written for Fast Fourier Transform (FFT), This was 
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oBsentially dona to conipare the F]?T results escpected with 
3-bits microprooessor with fixed -point arithmetic with 
rOvSults obtained with 36~bits 7044 coraputer with floating- 
point arithmetic so that a decision can be tahen about 
0 bits word length based on the results of comparison. 

4.4.1 Simulation of input data 

The input data will be normalized 5 bits plus sign, 
i.e. it xfill be between 011111 and 111111, where first 
bit is the sign bit and other 5 bits for magnitude and the 
binary point is assumed to lie to the left of leftmost 
magnitude bit. The input sequence (Real and imaginary) was 
generated to lie between 0 11111 and 1 11111 to simulate the 
input data. This in turn means generating: integers between 
+ 124 (0 mil 00) and - 124 (1 11111 00) for 8“bit machine. 
For 36-bits 7044 coraputer, the generated input sequence will 
be as shown in Pig. 4.3. 

4.4.2 7-Bits value of SIFF 

7-bits 257 values of SIrIB from 0 to rr/2 in step size 
of 2 ti;/ 1024 were calculated. Those values will vary between 
0000000 and 1111111. These values of SIEfE were read in to 
the simulated program. 

4.4.3 Simulation of 8— bits (7— bits for magnitude, on bit for 
sign) Machine 
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Sign bit 


^ 1 

1 

0 ,000000000 UOOOOOOOUOO . . . 00000 

1111100 

} 

i ^ ^ 

I 4 28 bits -> 

1 ( 

1 

i - Y "b i i5"S- 

* t 


Maximujm value of data + 124 


Sign lit 



! ' ' 

I 1 OOUOOOOOOOOOOOOOOOOO ... OOOuCOi 1111100 



28 "bits 

Minimum value of data "-124 



bit-s>^ 


Pig. 4.3 Data representation of S-lits maokine on 
7044 ( 36 -bits machine) 



¥lier>. two 7“bits (plus sign) numbers are multiplied, 
we jet a result of I 4 bits. Out of these 14 bits, MSB 
7 bits will be used in 8-bit imchine. In the simulated 
program., this was achieved by dividing the results of 
nultiplication, which are 14 “bits, by 2 ^ i.e, equivalent of 
Jiving 7 right shifts. 

4 . 4.4 Simulation of Fixed Point Ai’ithmetic 

As we iiave stated that wo aro goirg to use fixed point 
aritliraetic in the FFT program, the magnitudes of components 
calculated during different stages of computations should not 
exceed 7*-bits to prevent overflow. This can be achieved in 
the simulated program by ensuring that magnitude of members, 
calculated two at a time in the inner loop, does not exceed 
127 . If the absolute value of components is less than or 
equal to 127, the computations continue. However, when the 
absolute value of any component calculated exceeds 127, which 
is equivalent of having an over.’flow for 8-bits machine, the 
entire array is divided by two, scaling factor incremented by 
one and the members of that array are calculated again. 

4 . 4.5 Conversion Prom Integer to Beal 

In the fixed— point arithmetic, the results of computations 
will have integer values. The integer values are converted 
to their real equivalent as follows j 
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i) ]\!u] '':iply the integer velue by .5/2^, in general by 

"R— 1 

,5/2 , whore B is the nunber of bits used to represent 

raagnitude. This aessentially raeans assigning xfeights 

to the integer value represented in binaiT* The MSB 

t B"-l 

bit has a ireight ,5 and weight of LSB bit is .5/2 

ii) Multiply the above value by 2 to the power of scaling 
factor to get the final real value. 

4,4.6 Results of Simulation 

The simulated ITT FORTRAN program was run on IBM 7044 for 
1024 saiaplo points. Th- results of this program wore compared 
’.'Tith the results obtained by the FPT program using floating 
point arithmetic and 55— bits word length with the same input 
sequence data. The observations are given below ; 

i) vUt of 1024 components (coefficients), 896 became zero 
i.G. 87,5 percent of the result produced by the simulated 

program wore zero, 

ii) Other 128 components, which were not zeros, were not 
marginally out, but differ coiisiderably from the 
correct results. 

The above results are notj^all acceptable from accuracy 
point of view. Based on the results of simulation, it is 
concluded that 8-bits word length is not enough for FFT 
program employing fixed-point arithmetic. 
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4.5 16-Bits Simla tion on 70 H 

Since the 8 -bits word length did not prove enoijgh, 
the next step was to try for 16-bits word length, A PORTEAN 
program, which simlates 16-bits word length PPT program 
employing fixed point arithmetic, was written for 7044 computer. 
A data sequence as shown in Pig. 4.4 was generated to simulate 
the input data. 15 -bits 257 values of SIM were calculated and 
read into the program. Phis program employs the same tech- 
nique for simulation of 16-bits and fixed-point arithmetic as 
discussed for 3-bits simulation. The program was run on 7044 
computer for 1024 sample points. The results of the program 
were compared with the correct results and the observations 
are given below : 

i) 99.6 percent of components were within the range of 
,01 to .15 from their correct values 

ii) .4 percent of components differ by .16 to ,3 from their 
correct values 

Based on the above results, which are acceptable, 

16-bits word length is required for PPT program employing 
fixed -point arithmetic. As Veldh [ 6 ] has calculated for 
fixed -point Past Fourier Transform, the ratio of rms error 
to rms result with 16— bits word length will be of the order 
of 2 * lO""^. 
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Sign Bit 



I 

I 


Maximum Value of Data Generated ; + 31744 



Minimum Value of Data Generated i -31744 


Big. 4.4 Data Representation for 16 tits machine on 
7044 (36 bits machine) 



44 

4.6 Selection of 8~bits nicroprocessor Versus 16‘-bits 
Microprocessor 

Tl; c r 'Suits of 8-bits and 16-bits simulation have shown 
thot G-bits word length is not enough and as such 16-bits 
word length is required. There were two alternatives j 

i) either to select 16-bits microprocessor Tf© 9900 

and develop FPT software for it end thus trs;- to achieve 
tine constraint of 620 msecs for real-time application 
(problem specified) 

ii) or to go ahead with selected 8080 8-bits microprocessor 
(for which assembly language program Iie'S already been 
Coded) and use double-precision to achieve the accuracy. 
V/ith double-precision, it will not be possible to ■ 
achieve the time target for real-time application. 

The first alternative was not acceptable as TMS 9900 
microprocessor was not available and secondly there was no 
software support like cross-assembler, simulator available. 
Without software support, it would have resulted fn FPT 
program without being fully tested and debugged. This 
prograra could not have been used straightway in future when 
TMS 9900 would have been available. 

The second alternative was acceptable since there is 
software support available for testing and dehugging the 

ass enbly language program . And further , the tested program 
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can bo ru.. on available MDS-C'j kit or Mcro-78 c.'iaputer, both 
of which use 8080 microprocessor as their processor. This 
will result in proven PPT paclcage for available 8080 
microprocessor based machine* 

4*7 8080 Assembly Language Program with Double Precision 

Assembly language pro'^ran for PPT was rewritten to have 
double precision. In this program, the data word length is 
16 bits, 15 bite for riagnitude and one bit for sign. All 
subroutines which win bo called by PPT program were redesigned 
to work on data word length of 16 bits. IJvery care was taken 
in their coding so that they talce less amount of time for 
their execution. 

4.7.1 Testing of Assembly language Program 

The following sub-routines which will be called by PPT 
prograi*’. were taken first for jheir testing : 

i) Double Precision ffoltiply Routine (DPMI) 

This routine multiplies two signed 16 -bits numbers and 
produces signed 16— bits result. Negative result is 
repreoonted in 2’s complement. 

ii) Oomplement Routine ( . OOMPL) 

This routine converts the signed 16-bits numbers to 
their 2* s complement equivalent. 
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lii) Q'/orClovr Routine (OFL) 

ihi}.-> routine divides the 3.6— bits signed number by two 
i.c. nhiftG right by one bit. 

iv) Pov.'.^r of Tv/o Routine (L0rT2I) 

It finds the Log2l^" of number of sample (I) assuming 
N is povjr; r of two . 

v) Bit RovGrsa.l Routine (BITR) 

It converts the positive 16-bits number to its equiva- 
lent bit reversed number. This routine is required only 
in case of in-placc algorithm. 

vi) Unscrambling Routine (UNSCM) 

This routine orders the output sequence which will be in 
scrambled form in the co.so of in-place algorithm. 

The above assembly language routines were tested using 
the cro os-ass omblor and simua.ator for 8080 microprocessor 
available on 7044 computer with different combinations of 
data. Once these routines were proved correct, the different 
blocks of FFT program were proved for their correctness. 

The routines wore linked with the FPT program and the complete 
assembly language program was run on 7044 using cross- 
asBomblcr and simulator files. The program was tested 
for 8 sample points (N = 8) and the results obtained were 
compared with the actual results and were correct.. 



She further testing of the prograxa was done on Micro'- 
78 computer for different number of sample points, 

4.8 Micro -78 Implomentation 

Micro -78 computer uses 8080 micro processor as its 
processor and as such assembly language progrsin developed 
can be implemented straightway on this micro computer. 

For testing the assembly language program on Micro-78 
for different number of sample points, it is assumed that 
data is available on paper tape in hexadecimal signed magni- 
tude form. Eandom data (ZREAI, XIMG) were generated and 
punched on paper tape for N = 8, 128, 256, 512, 1024. As 
the actual data will be 5-bits plus sign, the data generated 
was between + 124 and - 124 in he::adecinal form, 

4.8.1 Sequence of Operation for Testing on Micro-78 

The main program calls ZTROM routine for clearing the 
lower bytes of XREAL and XIMAG. Lower bytes are made zeros 
since XRBALs and XIMAGs are two bytes long and the datas 
read will be one byte in length. The main program then 
calls READ routine, which reads the data from paper tape and 
then stores it in the ®Bs of XRBALs and XimOs. FFT program 
is then called which calculates the Fourier coefficients 
and stores them in XRBALs and XIMGs. Finally, the main 
program calls IRTFP routine. This routine converts the 
integer values of the coefficients into their reals 
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equivalent and prints them on to Tele-type. 

The program was tested for S' = 8, 128, 256, 512 and 
1024 sample points and all possible hugs were removed using 
the ODS-78 (on-line -debugging system) of Micro-78. ' The 
results were oompared with the results obtained with the 
simulated program run on 7044. The maximum difference 
between the results was of the order of .05 (for M = 1024). 
This error is caused due to integer to floating-point 
conversion, since the integers values of the coefficients 
obtained through Micro-78 and 7044 were some. 

4.8.2 Memory Requirement for PFT Implementation on Micro-78 

Memory requirement for FFT iraplementation consists 
of program area and data area as shown in Fig. 4.5. 

4. 8. 2.1 Program Area 

The program area consists of the following i 

i) The FFT program talces 1620 bytes of memory and has 
been assembled from 6000 to 11135. (octa-l) locations 
in memory of Micro -78, 

ii) The READ routine, which consists of IITPUT, ZEROM 

and READ routines, teJees 127 bytes and is assembled 
from 11150 to 11346 (octal) locations in neraory. 
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PROGRAM 


11135 

11150 


11346 

11350 


12002 

12010 

12057 

12060 


16057 


22057 



Pig. 4.5 Memory reciuirement for PFl implementation 
on Micro-78 
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iii) The INTIP routine, consisting of integer to floating- 
point conversion routine, Binary to BCD conversion 
routine (BNBCD) and printing routine (PRM!), takes 
305 bytes and has been assembled from 11350 to 12002 
(octal) locations. 

iv) The MAIH program, which calls the a,bove said routine, 
takes 40 bytes and is assembled from 12010 to 12057. 

Thus a total of 2090 bytes of memory area has been 
utilized for PPT and auxiliary prograi'as. 

4.8.2, 2 Data Area 

As the 16-bits word length is required, 21^ bytes for 
Real and 211 bytes for imaginary data points aro required for 
their storage (H is number of sample points). Thus a total of 
411 bytes of data area is required for FFT of N sample points, 

Dor IT = 1024, 4K bytes of data e.rea is required and for 
this va,luc of sample points, locations 12060 to 16057 have 
been reserved for real and locations 1S060 to 22057 for 
imaginary data points. 

4.9 Execution Tine 

The execution time for PET program for different number 
of sample points is given in Table 4.1. 

This time excludes the time taken for reading the data 
and the time for printing the results on the teletype. 
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Ho. of complex 
Sample points 


lime required for Time required for 

PPT progrem with PPT program with 

software nultiplica- hardware (8x8 hit) 

ti on routine multiplier 


128 


4 sec 500 millisec 2 sec 400 millisec 


256 10 sec 500 millisec 5 sec 600 millisec 

512 25 sec 200 millisec 12 sec 150 millisec 


1024 51 sec 600 millisec 28 sec 500 millisec 


Table 4.1 


Time required for the PPT program for 
different complex sample points with 
software multiplication routine and 
hardware (8x8 hit) multiplier. 
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4.10 Limitation of Program 

Tho complete program listing obtained tlirough Micro-78 
is given at Appendix ’P'. The program can be used for EFT of 
upto 1024 sample points. The limitation is because of SIME 
values stored in tho program, vhich are corresponding to 
N = 1024 points. The weights of ¥ (SHE and COSIHE) for 
lower number of sample points can be calculated from those 
values, but for number of samples greater thajx 1024, this is 
not possible. 


f 
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diAPTEE 5 

COEOLUSIOn .APTO PUTEilE ¥ORIC 

8080 microprocessor assembly Icnguase program developed 
and run on micro -78 computer wliicli uses 3080 microprocessor 
could not meet the real-time requ-irement of 1024 points FPT 
in 620 msecs. The time requirement could not be achieved 
because of 8080 being 8-bits microprocessor and the minimum 
requirement of 16 -bits word length (proved by simulation) and 
also due to nonavailability of multiplier chips to build 
16 * 16 bits comples multiplier unit. As j.t xms decided to 
use available 8-bits 8080 microprocessor due to the nonavai- 
lability of 16 bits microprocGSSor and/or its related soft- 
ware support, double precision was to be used in the EFT 
program for this microprocessor. The ^ early idea, of inter- 
facing of 8 3E 8 bits complex hardware multiplier with PPT 
program, with which we calcula,t8d 1000 msecs for 1024 points 
and estimated to a,chieve 620 msecs with faster version of 
8080, was to be abondonod. Double precision required 
16 * 16 bits complex hardware multiplier unit. The 
16 16 bits complex multiplier could not be built due to 

tho nonavailability of multiplier chips. As such software 
double precision multiply routine (to multiply 16 ± 16 bits 
number), which calls 8 8 bits loardware multiply unit was 

used in the FJ?T program. This resulted in 30 fold increase 
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in tile initial calculation of 1000 msecs for 1024 points FPT 
and gave us 30 sees for 1024 sample points. 

5.1 Conclusion 

The work done has resulted in the design of optamized 
PPT program (Chapter 3 ), which can bo coded in the assembly 
language of suitable 16~bits microprocessor like TD© 9900 to 
achieve the real-time environment in future. It has also 
proved that 8-bits word length is not enough and minimum of 
16-bits word length is required for PPT program using fixed- 
point arithmetic (result of simulation Chapter 4). And fur- 
ther, the 8080 assembly language program developed will be 
useful in the development of PPT assembly language program 
for 16-bits microprocessor or faster version of 8-bits micro- 
processor in future, finally, it has resulted in ^T soft- 
ware porkage for micro-78 and as such PPT upto 1024 sample 
points can be done on micro-computer instead on large computer, 

5.2 Idture Work 

The following related work to FPT is suggested to 
improve the perfornarce (Time wise) and flexibility of the 
developed FPT package ? 

i) In the developed PFT package, the input routine roads 

data from paper tapo. In actual practice, this data may 
bo required to bo read from tapo unit. It is suggested 
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that tape unit be interfaced with micro-TS computer 
and input routine be modified accordingly. 

ii ) SIHE values corresponding to sojnple points greater 

than 1024 be calculated and stored in ROMs. For example, 
if we calculate the SIKB value corresponding to 2048 
sample points and store in ROM, then plugging in this 
ROM in memory module will ma,ke the RI'T program applica- 
ble upto 2048 points. This' will remove the limitation 
of the developed program for its use upto 1024 points, 
because in the developed program SINE values correspond- 
ing to 1024 points have been stored. 

iii) 16 3£ 16 bits complex hardware multiplier be built when 
multiplier chips are available. This multiplier unit 
when interfaced with micro-78 will cut domi the execution 
time of the RET program considerably. With this and 
faster version of 8080, it may be possible to achieve 
real-time constraint. This will involve very little 
modifica.tion in the software. 

iv) l6-bits microprocessor like T® 9900 may be tried for 
Rx^T. This will eliminate the need of double precision 
in the program. If the software developed for this . 
microprocessor be interfaced with hardware multiplier 
unit, there will be great reduction in the execution 
tine of PPT (for real-time application). 
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v) Para-llcl processor organization instoad of sequential 
scheme may ho tried to exploit parallelism in the 
algorithm, 

vi) Higher "base algorithms (Base 4 > Base 8, Base 16) may he 
tried instoad of Base 2 algorithm to reduce computation 
time [1] * 
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Appendix 


FFT IN PLACE PROGRAM FOR UP TO 1024 SAMPLE 
N< SAMPLES) SHOULD BE POWER OF TWO 


006000 


ORG 

6000B 

006 0 0 0 

000 FFTIN: 

MOP 


006 0 01 

303 

JMP 

FFT 


164 

017 


257 VALUES OF SINE FROM 0 TO PIE/2 
IN STEP OF 2PIE/1024 STORED 
006004 000 SIN : DW 0B 

000 

006006 311 DW 31 IB 

000 

006010 1-56 2_2Z DW 622U 

# 0^00 3 . 


006012 133 
002 

006014 044 
003 

006016 355 

003 

006020 266 

004 

006022 177 

005 

006024 110 

006 

006026 021 
007 



011 

006036 063 
012 

006040 373 
012 

006042 304 

013 

006044 214 

014 

006046 124 

015 

006050 034 

016 

006052 344 
016 

006054 253 
017 

006056 163 
020 

006060 072 
021 

006062 001 
022 

006064 310 


DW 

1 133B 

DW 

1444B 

DW 

1755B 

DW 

2266B 

DW , 

2577B 

DW 

31 1 0B 

DW 

3421B 

DW 

3731B 

DW 

4242B 

DW 

4553B 

DW 

506 3B 

DW 

5373B 

DW 

5704B 

DW 

6214B 

DW 

6524E 

DW 

7034B 

DW 

7344B 

DW 

7653B 

DW 

10163B 

DW 

10472B 

DW 

1 1001B 

DW 

1 1310B 


022 


006066 

217 

023 

DW 

1 1617B 

006070 

125 

024 

DW 

12125B 

006072 

034 

025 

DW 

12434B 

006074 

342 

025 

DW 

12742B 

006 076 

250 

026 

D¥ 

13250B 

006100 

156 

027 

DW 

13556B 

006 1 02 

063 

030 

DW 

14063B 

006104 

371 

030 

DW 

14371B 

006 1 06 

276 

031 

DW 

14676B 

0061 10 

203 

032 

DW 

15203B 

006 1 12 

107 

033 

DW 

15507B 

006 1 1 4 

014 

034 

DW 

16fi4B 

006116 

320 

034 

DW 

16320B 

006120 

223 

035 

DW 

16623B 

006122 

127 

036 

DW 

17127B 

006124 

032 

037 

DW 

17432B 

006126 

335 

037 

DW 

17735B 

006130 

237 

040 

DW 

20237B 

006132 

142 

041 

DW 

20542B 

006134 

044 

042 

DW 

21044B 

006 1 36 

345 

042 

DW 

21345B 

006140 

247 

043 

DW 

21647B 

006142 

147 

044 

DW 

22147B 

006.144 

050 

045 

DW 

22450B 

006146 

350 

045 

DW . 

22750B 

006150 

250 

046 

DW 

23250B 

006152 

150 

DW 

23550B 



006154 

047 

04 7 

DW 

24047B 

006156 

050 

345 

DW 

24345B 

006160 

050 

244 

DW 

24644B 

006162 

051 

142 

DW 

25142B 

006164 

052 

037 

DW 

25437B 

006166 

053 

334 

DW 

25734B 

006170 

053 

231 

DW 

26231B 

006172 

054 

125 

DW 

26525B 

006174 

055 

021 

DW 

27021B 

006 1 76 

056 

314 

DW 

27314B 

006200 

056 

207 

DW 

276 07B 

006202 

057 

102 

DW 

30102B 

006204 

060 

374 

DW 

30374B 

006206 

060 

265 

DW 

30665B 

006210 

061 

its 6 

DW 

31156B 

006212 

062 

047 

DW 

31447B 

006214 

063 

337 

DW 

31737B 

006216 

063 

227 

DW 

32227B 

006220 

064 

116 

DW 

325 16B 

006222 

065 

004 

DW 

33004B 

006224 

066 

272 

DW 

33272B 

006226 

066 

160 

DW 

33560B 

006230 

067 

045 

DW 

34045B 

006232 

070 

331 

DW 

34331B 

006234 

070 

215 

DW' 

34615B 

006236 

071 

100 

DW 

351 00B 

006240 

072 

363 

, DW 

35363B 



006242 

072 

245 

006244 

073 

127 

006246 

074 

010 

006250 

075 

270 

006252 

075 

150 

006254 

076 

027 

006256 

077 

306 

006260 

077 

164 

006262 

100 

041 

006264 

101 

316 

006266 

101 

172 

006270 

102 

046 

006272 

103 

321 

006274 

103 

173 

006276 

104 

044 

006300 

105 

315 

006302 

105 

165 

006304 

106 

035 

006306 

107 

304 

006310 

107 

152 

006312 

110 

017 

006314 

1 1 1 
264 

006316 

11 1 
130 

006320 

1 12 
373 

006322 

112 

236 

006324 

1 13 
100 

006326 

1 14 
341 

006330 

114 

201 


D¥ 35645B 
DW 36127B 
DW 364l'0B 
DW 3667'0B 
DW 37150B 
DW 37427B 
DW 37706B 
DW' 40164B 
DW 4044 IB 
DW 40716B 
DW 41172B 
DW 41446B 
DW 4 172 IB 
DW 42173fe 
DW 42444B 
DW 42715B 
DW 43165B 
DW 43435B 
DW 43704B 
DW 44152B 
DW 444 17B 
DW 44664B 
DW 45130B 
DW 45373B 
DW 45636B 
DW 461 00B 
DW 46341B 
DW 4660 IB 



006332 

041 

116 

006334 

300 

116 

006336 

136 

117 

006340 

373 

117 

006342 

230 

120 

006344 

064 

121 

006346 

317 

121 

006350 

151 

122 

006352 

003 

123 

006354 

233 

123 

006356 

063 

124 

006360 

312 


124 

006362 

140 

125 

006364 

366 

125 

006366 

212 

126 

006370 

036 

127 

006372 

261 

127 

006374 

103 

130 

006376 

324 

130 

006400 

144 

131 

006402 

364 

131 

006404 

202 

132 

006406 

020 

133 

00641 0 

235 

133 

0064 1 2 

051 

134 

006414 

264 

134 

006416 

076 

135 

006420 

310 


D¥ 

47041B 

DW 

47300B 

D¥ 

47536B 

D¥ 

47773B 

D¥ 

50230B 

D¥ 

50464B 

D¥ 

507 17B 

D¥ 

51 15 IB 

D¥ 

51^038 

D¥ 

51633B 

D¥ 

52063B 

D¥ ' 

523 12B 

D¥ 

52540B 

DW 

52766B 

DW 

53212B 

DW 

53436B 

D¥ 

53661B 

DW 

54103B 

DW 

54324B 

DW 

54544B 

DW 

54764B 

DW 

55202B 

DW 

55420B 

DW 

55635B 

DW 

56051B 

DW 

56264B 

DW 

56476B 

DW 

567 10B 


135 

006422 120 

136 

006424 327 

136 

006426 136 

137 

006430 344 
137 

006432 150 
140 

006434 354 

140 

006436 157 

141 

006440 36 1 

141 

006442 162 

142 

006444 362 

142 

006446 161 

143 

006450 357 

143 

006452 154 

144 

1^6454 351 

144 

006456 144 

145 

006460 336 

145 

006462 127 

146 

006 4 6 4 32 0 

146 

006466 107 

147 

006 4 7 0 2 7 5 
147 

006472 062 
150 

006474 247 

150 

006476 032 

151 

006500 214 
151 

006 5 0 2 3 7 5 

151 

006504 156 

152 

006506 335 
152 

006 5 1 0 1 13 
15'3 


DW 

57120B 

DW 

57327B 

DW 

57536B 

DW 

57744B 

DW 

60150B 

DW 

60354B 

DW 

60557B 

DW 

60761B 

DW 

61 162B 

DW 

61362B 

DW 

61561B 

DW 

61757B 

DW 

62154B 

DW 

6235 IB 

DW 

62544B 

DW 

62736B 

DW 

63127B 

DW 

63320B 

DW 

63507B 

DW 

63675B 

DW 

64062B 

DW 

64247B 

DW 

64432B 

DW 

646 14B 

DW 

64775B 

DW 

65156B 

DW 

6 53 3 SB 

DW 

655 13B 



006 5 1 2 

270 

153 

006514 

044 

154 

006 5 1 6 

217 

154 

006520 

371 

154 

006522 

142 

155 

006524 

312 

155 

006526 

061 

156 

006530 

227 

156 

006532 

373 

156 

006534 

137 

•157 

006536 

302 

157 

00 6 5 4 0 

043 

160 

006542 

203 

160 

006544 

343 

160' 

00 6 5 46 

101 

161 

006550 

236 

161 

006552 

372 

161 

006554 

125 

162 

006556 

257 

162 

006 56 0 

010 

163 

006562 

137 

163 

006564 

266 

163 

006 5 6 6 

013 

164 

006 5 7 0 

140 

164 

006572 

263 

164 

006574 

005 

165 

006576 

126 

165 

006600 

246 


DW 

65670B 

D¥ 

66044B 

DW 

66217B 

DW 

66371B 

DW 

66542B 

DW 

667 12B 

DW 

6706 IB 

DW 

67227B 

DW 

•67373B 

DW , 

67537B 

DW 

67702B 

DW 

70043B 

DW 

70203B 

DW 

70343B 

DW 

70501B 

DW 

70636B 

DW 

70772B 

DW 

7 1 125B 

DW 

7 1257B 

DW 

714 10B 

DW 

71537B 

DW 

7 1666B 

DW 

720 13B 

DW 

72140B 

DW 

72263B 

DW 

72405B 

DW 

72526B 

DW 

72646B 



165 

006602 364 

165 

006604 102 

166 

006606 216 
166 

006610 331 
166 

006612 043 
167 

006614 154 
167 

006616 264 
167 

006620 373 
167 

006622 100 
170 

006624 205 
170 

006626 310 

170 

006630 012 

171 

006632 112 
* 171 
006634 212 

171 

006636 311 
17 1 

006640 006 

172 

006642 102 
172 

006644 175 
172 

006646 267 
172 

006650 357 

172 

006652 047 

173 

006654 135 
173 

006656 222 
173 

006660 306 
173 

006 6 6 2 3 71 

173 

006 6 6 4 0 5 2 

174 

006666 132 
174 

006 6 7 0 21 1 


D¥ 

72764B 

D¥ 

73102B 

D¥ 

73216B 

D¥ 

73331B 

D¥ 

73443B 

D¥ 

73554B 

D¥ 

73664B 

D¥ 

73773B 

D¥ 

74100B 

D¥ 

74205B 

D¥ 

743 10B 

D¥ 

744 12B 

D¥ 

745 12B 

D¥ 

746 12B 

D¥ 

7471 IB 

D¥ 

75006B 

D¥ 

75102B 

D¥ 

75175B 

D¥ 

75267B 

D¥ 

75357B 

D¥ 

75447B 

D¥ 

75535B 

D¥ 

75622B 

D¥ 

75706B 

D¥ 

7577 IB 

D¥ 

76052B 

D¥ 

76132B 

D¥ 

7621 IB 



174 

306672 267 
174 

306674 344 

174 

306676 017 

175 

ZH36700 072 
175 

306702 143 
175 

006704 212 
175 

006706 261 
175 

006710 326 
175 

006712 373 

175 

006714 036 

176 

006716 077 
176 

006720 140 
176 

006722 177 
176 

006724 235 
176 

006726 272 
176 

006730 326 
176 

006732 360 

176 

006734 012 

177 

006736 042 
177 

0067 4 0 0 7 0 
1.77 

006742 116 
177 

006744 142 
177 

006746 165 
177 

006 7 5 0 2 0 5 
177 

006752 230 
177 

P6754 247 
177 

006 7 5 6 2 6 5 
177 

006 7 6 0 3 02 


DW 

76267B 

DW 

76344B 

DW 

76417B 

DW 

76472B 

DW 

76543B 

DW 

766 12B 

DW 

76661B 

DW 

76726B 

DW 

76773B 

DW 

77036B 

DW 

77077B 

DW 

77140B 

DW 

77177B 

DW 

77235B 

DW 

77272B 

DW 

77326B 

DW 

77360B 

DW 

77412B 

DW 

• 77442B 

DW 

77470B 

DW 

77516B 

DW 

77542B 

DW 

77565B 

DW 

77605B 

DW 

77630B 

DW 

77647B 

DW 

77665B 

DW 

77702B 



006762 

177 

316 


DW 

77716B 

006 764 

177 

331 


D¥ 

7773TB 

mei66 

177 

342 


DW 

77742B 

006770 

177 

352 


DW 

77752B 

006 7 7 2 

177 

361 


DW 

77761B 

meiiA 

177 

366 


DW 

77766B 

006776 

177 

372 


DW 

77772B 

00 7 0 00 

177 

376 


DW 

77776B 

007 0 0 2 

177 

377 


DW 

77777B 

007004 

177 

377 


DW 

77777B 

007 006 

177 

000 

NI2 s 

DW 

0B 

00 7 0 1 0 

000 

000 

I : 

DW 

0B 

007012 

000 

000 

N : 

DW 

0B 

a570t4 

000 

000 

lA t 

DW ' 

0B 

00 7 016 

000 

000 

LL ; 

DB 

0B 

007017 

000 

HH : 

DB 

0B 

007020 

000 

CT512 : 

DW 

1000B 

007022 

002 

000 

S1024: 

DW 

0B 

007024 

000 

000 

C1024J 

DW 

2000B 

007026 

004 

000 

SI512: 

DW 

0B 

mimm 

000 

000 

SN512 ; 

DW - 

0B 

00 7 0 3 2 

000 

000 

C512 : 

DW 

177000B 

376 

eeae 0 

007034 001 

IB : 

DW 

IB 

007036 

000 

000 

IC ; 

DW 

0B 

007040 

000 

000 

ID ; 

DW 

0B 

007 0 4 2 

000 

000 

I2D : 

DW 

0B 

007 0 44 

000 

000 

K : 

DW 

0B 

007046 . 

000 

000 

SINE : 

DW 

0B 


000 





OUTER LOOP PARAMETER 
; NUMBER OF SAMPLE POINTS 
INITIALLY WILL HAVE N/2 
STORES L0G2N 



0B 


007050 000 
000 

007052 000 
007 0 5 3 0 0 0 
007054 000 
007055 000 
007056 000 
000 

00 7 0 6 0 0 0 0 
000 

007062 000 
000 

007064 000 
000 

007066 000 
000 

007 0 7 0 0 0 0 
000 

007072 000 
000 

007074 000 
000 

007076 000 
000 

007100 000 
000 

007102 000 
000 

mi\0ii 000 
,000 
007106 000 
000 

007110 000 
000 

007112 000 
000 

007114 000 
000 

0071 16 000 
007117 000 
007120 000 
007121 257 


007 1 2 2 06 2 

053 
016 

007 1 2 5 0 6 2 

054 
0L6 

007130 174 
007 1 3 1 3 46 
200 

007 1 3 3 3 1 2 
156 


COSIN; DW 

COUNT: DB 0B 

FLAG 1 ; DB 0B 

FLAG2 : DB 0B 

FLAG ; DB 0B 

EL : DW 0B 

EH ; DW 0B 

DL : DW 0B 

DH ; DW 0B 

MIR ; DW 0B 

Mil : DW 0B 

NMR : DW . 0B 

NMI : DW 0B 

XYRMI ; DW 0B 

XYIMI : DW 0B 

AREAL: DW 0B 

AIMAG: DW 0B 

REALl: DW 0B 

BREAL : DW . 0B 

IMAG 1 : DW 0B 

BIMAG: DW 0B 

FLAG4 : DB 0B 

BFLAG ; DB 0B 

LL8 : DB 0B 

DPMUL; XRA A 

DPMUL multiplies 16*16 BITS SIGNED NUMBER 

negative number in twos COMPLEMENT 

WHEN CALLED TWO NUMBERS IN CD^E) AND (H>L> REG* 
RESULT IN (H^L) REG. 

STA FLAGl 


STA FLAG2 


MOV A>H 
ANI 200B 

JZ DCHEC 




016 




007 1 36 

076 


MV I 

A^ IB 


001 




007 1 40 

062 


STA 

FLAGl 


053 





016 




007143 

257 


XRA 

A 

007144 

175 


MOV 

A^L 

037145 

057 


CMA 


037146 

306 


ADI 

IB 


001 




037150 

157 • 


MOV 

L^A 

g®7151 

174 


MOV 

AjH 

007 1 52 

057 


CMA 


007153 

316 


AC I 

0B 


000 




007 1 55 

147 


MOV 

H^ A 

007 1 56 

172 

DCHECs 

MOV 

A-.D 

007 1 57 

346 


ANI 

200B 


200 




e®7161 

312 


JZ 

PROdD 


204 





016 




007 164 

076 


MVI 

A^ IB 


001 




007166 

062 


STA 

FLAG2 


054 





016 




e®7171 

257 


XRA 

A 

007 1 72 

173 


MOV 

A*E 

007173 

057 


CMA 


007 1 74 

306 


ADI 

IB 


001 




007176 

137 


MOV 

E^A 

007 1 77 

172 


MOV 

A^D 

037 2 0 0 

057 


CMA 


007201 

316 


AC I 

0B 


000 




00 7 2 0 3 

127 


MOV 

D^A 

00 7 2 0 4 

072 

PROCD: 

LDA 

FLAG 1 


053 





016 




007207 

107 


MOV 

B^A 

007210 

072 


LDA 

FLAG2 


054 





016 




007213 

250 


XRA 

B 

007 2 1 4 

062 


STA 

FLAG 


055 





016 




007217 

325 


PUSH 

D 

007220 

345 


PUSH 

H 

00722 1 

1 14 


MOV 

C^H 

007222 

175 


MOV 

AjL 

007223 

315 


CALL 

MUL 


345 





016 






007226 

016 

007231 171 
007232 315 
345 
016 

007235 042 
060 
016 

007240 341 
007241 321 
007242 175 
007243 132 
007244 315 
345 
016 

007247 042 
062 
0 16 

007252 171 
007253 315 
345 
016 

007256 042 
064 
016 

007261 052 
060 
016 

007264 353 
007265 052 
062 
016 

007270 031 
00 7 2 7 1 3 5 3 
007272 052 
056 
016 

007 2 7 5 1 54 
®37276 ' 046 
000 

2837300 031 
037301 175 
007 3 0 2 3 5 3 
037303 132 
037304 026 
000 

037306 052 

064 

016 

M731 1 031 
007312 051 
007313 027 
007314 175 


SHLD 

EL 

MOV 

A^C 

CALL 

MUL 

SHLD 

EH 

POP 

H 

POP 

D 

MOV 

A#L 

MOV 

E>D 

CALL 

MUL 

SHLD 

DL 

MOV 

AaC 

CALL 

MUL 

SHLD 

DH 

LHLD 

EH 

XCHG 

LHLD 

DL 

DAD 

D 

XCHG 

LHLD 

EL 

MOV 

L^H 

MV I 

H>0B 

DAD 

D 

MOV 

A/L 

XCHG 

MOV 

E/D 

MV I 

Dv 0B 

LHLD 

DH 


DAD 

D 

DAD 

H 


RAL 

MOV 



007315 

006 


MVl 

B^0B 


000 




007 3 1 7 

210 


ADC 

B 

007 3 2 0 

157 


MOV 

L>A 

007321 

072 


LDA 

FLAG 


055 





016 




00 7 3 24 

376 


CPI 

0B 


000 




037326 

312 


JZ 

- FINIl 


344 





016 




007331 

257 

# 

XRA 

A 

007332 

175 


MOV 

A^L 

007333 

057 


CMA 


207334 

306 


ADI 

IB 


001 




037 3 3 6 

157 


MOV 

L^A 

007337 

174 


MOV 

A/H 

00 7 3 40 

057 


CMA 


00 7 341 

316 


AC I 

0B 


000 




00 7 3 4 3 

147 


MOV 

H>A 

007344 

311 

FINIl : 

RET 


007 3 4 5 

041 

MUL : 

LXI 

H>0B 


000 





000 




£07350 

006 

010, 


MVI 

10B 

007352 

026 


MV I 

D^0B 


000 




007354 

051 

MULLP; 

DAD 

H 

007355 

027 


RAL 


007356 

322 


JNC 

DEC . 


364 





016 




007361 

031 


DAD 

D 

007362 

316 


AC I 

0B 


000 




007364 

005 

DEC : 

DCR 

B 

007 3 6 5 

302 


JNZ 

MULLP 


354 





016 




00 7 3 7 0 

311 


R:ET 



COMPL FINDS TWOS COMPLEMENT OF 16 BITS NUMBER 
WHEN CALLED NUMBER IN <H>L> REG. 

RESULT IN <HaL) REG. 


00 7 3 7 7 


ORG 

7377B 

007377 

000 COMPL: 

NOP 


007400 

257 

XRA 

A 

007401 

175 

MOV 

A^L 

007402 

057 

CMA 


007403 

306 

ADI 

01B 

,♦ ■ ' V 

001 



00 7 4 0 5 

157 

MOV 

L'A 

007406 

174 

MOV , 

A^H 

007 4 0 7 

057 

CMA 



00 7 4 1 0 

316 

AC I 

0B 


000 



C®7412 

147 

MOV 

H«A 

00 7 4 1 3 

311 

RET 



OFL DIVIDES 16 BITS SIGNED NUMBER BY 2 
WHEN CALLED NUMBER IN (H^L) REG. 




RESULT 

IN (Hj 

L) REG 

007 4 1 4 

000 

OFL : 

NOP 


00 7 4 1 5 

076 


MV I 

A^0B 

- 

000 




007417 

062 


STA 

FLAG 4 


116 





016 




IZ07422 

174 

- ' 

MOV 

A,H 

e®7423 

346 


AN I 

200B 


200 




00 7 4 2 5 

312 


JZ 

DVIDE 


040 





017 




00 7 4 3 0 

076 


MVI 

A^ IB 


00 1 




007 4 32 

062 


STA 

FLAG4 


1 16 





016 




007435 

315 


CALL 

GOMPL 


377 





016 



■Sv, ■■■■ 

007440 

000 

DVIDE:' 

NOP 


00 7 441 

257 


XRA 

A 

007442 

174 


MOV 

A>H 

007443 

037 


RAR 


00 7 4 44 

147 


MOV 

H^A 

007 4 4 5 

175 


MOV 

A^L 

00 7 4 4 6 

037 


RAR 


007447 

157 


MOV 

LjA 

00 7 4 5 0 

072 


LDA 

FLAG 4 


1 16 





016 




007453 

376 


CPI 

IB 


001 




007 4 5 5 

372 


JM 

OVER 


063 





017 




007460 

315 


CALL 

COMPL 


377 





016 




007463 

31 1 

OVER : 

RET 



B I TR CALCULATES BIT REVERSED NUMBER OF 16 BITS NUMBER 
WHEN CALLED NUMBER INCHED REG. 

RESULT IN <H>L) REG. 


007464 

000 BITR ! 

NOP 


007465 

353 

XCHG 


007466 

072 

LDA 

BFLAG 


1 17 




016 



00747 1 

376 

CPI 

IB 




001 




027473 

372 


jm 

LLLE8 


137 

017* 




007476 

257 

LLGT8: 

XRA 

A 

miaii 

173 


MOV 

A^E 

007500 

036 


MVI 

E^ 10B 


010 




007502 

041 


LXI 

H^00B 


000 





000 




027505 

051 

BRLP3 t 

DAD 

H 

007 5 06 

037 


RAR 


007 5 0 7 

322 


JNC 

BRLP4 


113 





017 




027512 

043 


INX 

H 

007 5 1 3 

035 

BRLP4 : 

DCR 

E 

007 5 1 4 

302 


JNZ 

BRLP3 


105 





017 




00 7 5 1 7 

072 

120 


LDA 

LL8 


016 




007 5 2 2 

137 


MOV 

E>A 

007 5 2 3 

172 


MOV 

A^ D 

007524 

051 

BRLP5 s 

DAD 

■M'/ 

e®7525 

037 


RAR 


007526 

322 


JNC 

BRLP6 


132 





017 




02 7 5 31 

043 


INX 

H 

007532 

035 

BRLP6 : 

DCR 

E 

007533 

302 


JNZ 

BRLP5 


124 





017 




007536 

311 


RET 


007537 

000 

LLLE8S 

NOP 


007540 

257 


XRA 

A 

007541 

041 


LXI 

H^LL 


016 





016 




007544 

173 


MOV 

A^E 

007545 

136 


MOV 

E>M 

007546 

041 


LXI 

H#0B 


000 





000 




007551 

051 

BRLPl : 

DAD 

H 

007552 

037 


RAR 


007 5 5 3 

322 


JNC 

BRLP2 


157 





017 




007 5 5 6 

043 


INX 

H 

00 7 5 5 7 

035 

BRLP2 : 

DCR 

E 

007 5 6 0 

302 


JNZ 

BRLPl 


151 







017 




M7563 

31 1 

RET 




FFT PROGRAM 

Starts from 

HERE 

007564 

052 FFT s 
012 

016 

LHLD 

N 


007 5 6 7 

315 

■mt io3 

022 

CALL 



00 7 5 7 2 

m 

257 

XRA 

A 


mis? 3 

174 

MOV 

A^H , 


007574 

037 

RAR 



e®7575 

147 

MOV 

H#A 


misie 

175 

MOV 

A, L 


misn 

037 

RAR 



00 7 6 0 0 

157 

MOV 

*L>A 


007 6 01 

042 

SHLD 

lA 



014 





016 




00 7 6 04 

042 

006 

SHLD 

NI2 



016 




007607 

076 

MVI 

A^ 1 0B 



010 




00761 i 

270 

CMP 

B 



372 

JM 

GTS ITO 

CHECK L0G2N >80R <8 


224 





017 




00 7 6 1 5 

257 

XRA 

A 


007616 

062 

STA 

BFLAG 



117 





016 




00 7 6 21 

303 

JMP 

LTGT8 



237 





017 




037624 

076 GT8 : 
001 

MVI 

A^ IB 


007626 

062 

STA 

BFLAG 



117 





016 




00 7 6 3 1 

170 

MOV 

AyB 


00 7 6 3 2 

326 

SUI 

10B 



010 


' 


00 7 6 3 4 

062 

STA 

LL8 



120 





016 




007637 

170 LTGT8: 

MOV 

A^B 


007640 

062 

016 * 

016 

STA 

LL 


007643 

326 

SUI 

02B 



002 




007645 

107 

MOV 

B>A 


007646 

257 

XRA 

A 


007647 

076 

MVI 

A^01B 




007 6 51 

001 

037 

RAR 


00 7 6 52 

037 

LOOP : RAR 


007 6 5 3 

005 

OCR 

B 

00 7 6 54 

302 

JNZ 

LOOP 

007 6 5 7 

252 

017 

062 

STA 

HH 

00 7 6 6 2 

0 1 7 
016 
052 

LHLD 

Cl 024 

miees 

024 

016 

021 

Lxr 

D^ SIN 


004 

014 

031 

DAD 

D 

0076 7 1 

042 

SHLD 

S I 024 

mibiti 

022 

016 

052 

LHLD 

C 5 1 2 

mien 

032 

016 

353 

XCHG 


mnm 

041 

LXI 

H/SIN 

mnez 

004 

014 

031 

DAD 

D 

007704' 

042 

SHLD 

SI512 

mim 

026 

016 

052 

LHLD 

CT512 

mn 12 

020 

016 

021 

LXI 

D^SIN 

mil 1 5 

004 

014 

031 

DAD 

D 

007716 

042 

SHLD 

SN512 

007721 

030 

016 

257 

XRA 

A 

007722 

062 

STA 

I 

00 7 7 2 5 

010 

016 

041 

" LXI 

H^01B 

mn30 

001 

000 

042 

' SHLD 

IB 

007733 

034 

016 

04 1 

OUTER LOOP STARTS 
OUTER: LX;i H*0B 


000 

000 





SHLD IC 


007736 042 
036 
016 

mn4l 042 SHLD K 

044 
016 

007744 052 LHLD lA 

014 
016 

m77A7 042 SHLD ID 

040 
016 

007752 051 DAD H 

007753 042 SHLD I2D 

042 

007756 052 INNERs LHLD lA ;INNER LOOP STARTS 

014 
016 

007761 353 XCHG 

007762 052 LHLD IC 

036 
016 

007765 000 ICIA : NOP 

POWER OF W CALCULATED IN HERE 


007 766 

257 


XRA 

A 

007767 

172 


MOV 

A^D 

007 7 7 0 

037 


RAR 


007771 

127 


MOV 

D>A 

007 7 7 2 

173 


MOV , 

A>E 

007773 

037 


RAR 


007774 

137 


MOV 

E#A 

007775 

332 


JC 

OICIA 


011 





020 




010000 

174 


MOV 

A>H 

010001 

037 


RAR 


010002 

147 


MOV 

H> A 

010003 

175 


MOV 

A>L 

010004 

037 


RAR 


010005 

157 


MOV 

L^A 

010006 

303 


JMP 

ICIA 


365 

1 




017 




01001 1 

000 

OICIA: 

NOP 


0100 12 

315 


CALL 

BITR 


064 





017 




010015 

257 


XRA 

A 

010016 

072 


LDA 

HH 


017 





016 




010021 

037 

I CLP : 

RAR 


010022 

376 


CPI 

0B . 


000 




010024 

312 


JZ 

FINIS 




033 




020 



010027 

051 

DAD 

H • 

010030 

303 

JMP 

ICLP 


021 




020 



010033 

000 

FINIS: NOP 




CHECK POWER 

OF W .> OR <= 256 

010034 

000 

GL256: NOP 


010035 

174 

MOV. 

A^H 

010036 

376 

CPI 

IB 


00 1 



010040 

332 

JC 

LT257 

. 

133 




020 



010043 

376 

CPI 

IB 


001 



010045 

302 

JNZ 

GT256 


056 




020 



010050 

175 

ACHEC : MOV 

A>L 

010051 

376 

CPI 

IB 


001 



010053 

332 

JC 

' LT257 


133 




020 





VALUES OF SINE AND COSINE CALCULATED FOR W>256 

010056 

000 

GT256 : NOP 


010057 

051 

DAD 

H 

010060 

353 

XCHG 


010061 

052 

LHLD 

S1024 


022 




016 



010064 

257 

XRA 

A 

010065 

1 75 

MOV 

A^L 

010066 

223 

SUB 

E 

010067 

157 

MOV 

L^ A 

010070 

174 

MOV 

A^H 

010071 

232 

SBB 

D 

010072 

147 

MOV 

HjA 

010073 

176 

MOV 

A^M 

010074 

001 

LXI 

B^SINE 


046 




016 



010077 

002 

STAX 

B 

010 1 00 

043 

INX 

H 

010101 

003 

INX 

B 

010102 

176 

MOV 

A/M 

010103 

002 

STAX 

B 

010104 

052 

LHLD 

SI512 


026 




016 



010107 

031 

DAD 

D 

0101 10 

001 

LXI 

B^COS.IN 


050 





016 


0101 13 

257 

XRA 

A 

0101 14 

176 

MOV 

A^M 

0101 15 

057 

CMA 


0101 16 

306 

ADI 

IB 


001 



010120 

002 

STAX 

B 

010121 

043 

INX 

H 

010122 

003 

INX 

B 

010123 

176 

MOV 

A#M 

010124 

057 

CMA 


010125 

316 

AC I 

0B 


000 



010127 

002 

STAX 

B 

010130 

303 

JMP 

INI 


176 

020 


VALUES OF SIMF ANB COSINE CALCULATED FOR POWER OF ¥=<256 


010133 

000 

LT257 s NOP 


010134 

051 

DAD 

H 

010135 

021 

LX I 

D^SIN 


004 




014 



010140 

353 

XCHG 


010141 

031 

DAD 

D 

010142 

001 

LXI 

B^SINE 


046 




016 



010145 

176 

MOV 

A^M 

01 0 1 46 

002 

STAX 

B 

010147 

043 

INX 

H 

010150 

003 

INX 

B 

010151 

176 

MOV 

A>M 

010152 

002 

STAX 

B 

010153 

052 

LHLD 

SN5 12 


030 




016 



010156 

257 

XRA 

A 

01 0 1 57 

175 

MOV 

A/L 

010160 

223 

SUB 

E 

010161 

157 

MOV 

L^A 

010162 

174 

MOV 

A^H 

010163 

232 

SBB 

D 

010164 

147 

MOV 

H>A 

010165 

176 

MOV 

A^M 

010166 

001 

LXI 

B^COSIN 


050 




016 



010171 

002 

STAX 

B 

010172 

043 

■ INX 

H 

010173 

003 

INX 

B 

010174 

176 

MOV 

A^M 

010175 

002 

STAX 

B , 

010176 

000 

INI : NOP 


010177 

052 

LHLD 

lA 




014 




016 



010202 

051 

DAD 

H 

010203 

353 

XCHG 


010204 

041 

LXI 

H^XREAL 


060 




024 



010207 

031 

DAD 

D 

010210 

042 

SHLD 

MIR 


066 




016 



010213 

041 

LXI 

H# XI MAG 


060 




034 



010216 

031 

DAD 

D 

010217 

042 

SHLD 

MI I 


070 




016 



010222 

052 

LHLD 

ic 


036 




016 



010225 

051 

DAD 

H 

010226 

115 

MOV 

C^L 

010227 

104 

MOV 

B/H 



INNER MOST LOOP STARTS 

010230 

052 

IMOST: LHLD 

MIR IBREAL CALCULATED 


066 




016 



010233 

01 1 

DAD 

B 

010234 

305 

PUSH 

B 

010235 

136 

MOV 

E^M 

010236 

043 

INX 

H 

010237 

126 

MOV 

D>M 

010240 

353 

XCHG 


010241 

315 

CALL 

OFL ;to prevent QEKFLOW 


014 




017 





mGE 001 


010244 


ORG 

10244B 

010244 

042 

076 

016 

SHLD 

XYRMI 

010247 

353 

XCHG 


010250 

052 

050 

016 

LHLD 

COSIN 

01 0253 

315 

121 

016 

CALL 

DPMUL 

010256 

042 

106 

016 

SHLD 

REALl 

01026 1 

301 

POP 

B 

010262 

052 

070 

016 

LHLD 

Mil 

010265 

011 

DAD 

B 

010266 

136 

MOV 

E>M 

010267 

043 

INX 

H 

010270 

126 

MOV 

D^M 

01027 1 

353 

XCHG 


010272 

315 

014 

017 

CALL 

OFL 

010275 

042 

100 

016 

SHLD 

XYIMI 

010300 

353 

XCHG 


010301 

052 

046 

016 

LHLD. 

SINE 

010304 

305 

PUSH 

B 

010305 

315 

121 

016 

CALL 

DPMUL 

f 

0103 10 

353 

XCHG 


0103 1 1 

052 

106 

016 

LHLD 

REAL 1 

0103 14 

03 1 

DAD 

D 

0103 15 

042 

1 10 
016 

SHLD BREAL 

BIMAG CALCULATED 

010320 

052 

100 

016 

LHLD 

XYIMI 

010323 

353 

XCHG 


010324 

052 

050 

016 

LHLD 

COSIN 

010327 

315 

121 

016 

CALL 

DPMUL 

010332 

042 
112 ' 
016 

SHLD 

I MAG 1 


010335 052 


LHLP XYRMI 




016 

010340 

353 

010341 

052 


046 


016 

010344 

315 


121 


016 

010347 

353 

010350 

052 


1 12 


016 

010353 

353 

010354 

315 


377 


016 

010357 

353 

010360 

031 

010361 

042 


1 14 


016 

010364 

301 

010365 

041 


060 


024 

010370 

01 1 

010371 

136 

010372 

043 

010373 

126 

010374 

353 

010375 

315 


014 


017 


042 


102 



010403 

041 


060 


034 

010406 

01 1 

010407 

136 

010410 

043 

01041 1 

126 

010412 

353 

01 04 1 3 

315 


014 


017 

010416 

042 


10A 

016 

SUM OF 

010421 052 
102 
016 

010424 353 
010425 052 
110 
016 

010430 031 


XCHG 


LHLD 

SINE 

CALL 

DPMUL 

XCHG 

LHLD 

IMAGl 

XCHG 

CALL 

COMPL 

XCHG 


DAD 

D 

SHLD 

BIMAG 

POP 

B 

LXI 

H..XREAL 

DAD 

B 

MOV 

E^M , 

I NX 

H' 

MOV 

XCHG 

D^M 

CALL 

OFL 

SHLD 

AREAL 

LXI 

H^XIMAG 

DAD 

B 

MOV 

; E^M 

I NX 

H 

MOV 

XCHG 

D^M 

CALL 

OFL 

SHLD 

AIMAG 

AREAL 

AND BREALCREAL 

LHLD 

AREAL 

XCHG 

LHLD 

BREAL 

DAD 

D 


VALUE OF COMPONENT) 


CALCULATE 


010431 

353 

XCHG 


010432 

041 

LXI 

H^XREAL 


060 




024 



010435 

011 

DAD 

B 

010436 

163 

MOV 

M^E 

010437 

043 

I NX 

H 

010440 

162 

MOV 

M#D 



AREAL -BREAL 'i 

CALCULATED 

010441 

052 

LHLD 

AREAL 


102 




016 



010444 

353 

XCHG 


010445 

052 

LHLD 

BREAL 


110 




016 



010450 

315 

CALL 

COMPL 


377 




016 



010453 

353 

XCHG 


010454 

031 

DAD 

D 

010455 

353 

XCHG 


010456 

052 

LHLD 

MIR 


066 




016 



010461 

01 1 

DAD 

B 

010462 

163 

MOV 

Mj.E 

010463 

043 

INX 

H 

010464 

162 

MOV 

MjD 

010465 

052 

LHLD 

.AIMAG 


104 




016 



010470 

353 

XCHG 


010471 

052 

LHLD 

BIMAG 


1 14 




016 



010474 


DAD 

D 

010475 

353 

XCHG 


010476 

041 

LXI 

H^XIMAG 


060 




034 



010501 

01 1 

DAD 

B 

01 0502 

163 

MOV 

M>E 

010503 

043 

INX 

H 

010504 

162 

. MOV 

M^D 

010505 

052 

LHLD 

AIMAG 


104 




016 



010510 

353 

XCHG 


01051 1 

052 

LHLD 

BIMAG 


1 14 




016 



010514 

315 

CALL 

COMPL 


377 




016 



010517 

353 

XCHG 

' - 

010520 

031 

DAD 

D 

010521 

353 

XCHG 


010522 

052 

LHLD 

Mil- 


070 




;aimag+bimag 


;aimag-bimag 


CALCULATED 


CALCULATED 



016 

010525 01 1 
010526 163 
010527 043 
010530 162 
010531 003 
010532 003 
'010533 171 
010534 041 
042 
016 

010537 276 
010540 302 
230 
020 

010543 043 
010544 170 
010545 276 
010546 302 
230 
020 

010551 052 
014 
016 

010554 051 
010555 353 
010556 052 
036 
016 

010561 031 
010562 042 
036 
016 

010565 052 
040 
. 016 
010570 031 
010571 042 
040 
016 

010574 051 
010575 042 
042 
016 

010600 052 
044 
016 

010603 043 
010604 042 
044 
0 16 

010607 353 
010610 041 
034 
016 

010613 173 
010614 276 
010615 302 
356 
017 . 


DAD 

B 

MOV 

M^E 

I NX 

H 

MOV 

M>D 

I NX 

B 

I NX 

B 

MOV 

A^C 

LXI 

H#I2D 


CMP 

M 

JNZ 

IMOST 

I NX 
MOV 

H 

A^B 

CMP 

M 

JNZ 

IMOST 

LHLD 

lA . 

DAD 

H 

XCHG 

LHLD 

IC 

DAD 

D 

SHLD 

IC 

LHLD 

ID 

DAD 

D 

SHLD 

ID 

DAD 

H 

SHLD 

I2D 

LHLD 

K 

I NX 

H 

SHLD 

K 

XCHG 

LXI 

H> IB 

MOV 

A^E 

CMP 

M 

JNZ 

INNE 



010620 043 

INX 

H 

010621 172 

MOV 

A>D 

010622 276 

CMP 

M 

010623 302 

JNZ 

INNER 

356 



017 



010626 257 

XRA 

A 

010627 052 

LHLD 

IB 

034 



016 



010632 175 

MOV 

A>L 

010633 027 

RAL 


010634 157 

MOV 

L^A 

010635 174 

MOV 

A^H 

010636 027 

,HAI- 


010637 147 

MOV 

H^A 

010640 042 

SHLD 

IB 

034 



016 



010643 257 

XRA 

A 

010644 052 

LHLD 

lA 

014 



016 



010647 174 

MOV 

A^H 

010650 037 

RAR 


010651 147 

MOV 

H^A 

010652 175 

MOV 

A#L 

010653 037 

RAR 


010654 157 

MOV 

L>A 

010655 042 

SHLD 

lA 

014 



015 



010660 072 

LDA 

I 

010 



016 




INR 

A 


STA 

I 

010 



016 



010667 041 

LX I 

H^LL 

016 



016 



010672 276 

GMP 

M 

010673 302 

JNZ 

OUTER 

333 



017 




UNSCRAMBLING 

OF OUTPUT 

010676 001 

LXI 

B » 0B 

000 



000 

- 


010701 151 

XLPl ; MOV 

LjC 

010702 140 

MOV 

H^B 

010703- 305 

PUSH 

B 

010704 315 

CALL 

BITR 

064 



017 



010707 174 

MOV 

A^H 

0107 10 270 

CMP 

B 

010711 372 

JM 

LPCON 

, 054 





CHECA 


022 

010714 312 
041 

» 022 
010717 000 EXCHG 
010720 051 
010721 042 
134 
022 

010724 353 
010725 151 
010726 140 
010727 051 
010730 115 
010.731 104 ■ 

010732 041 
060 
024 

010735 031 
010736 136 
010737 043 
010740 126 
010741 053 
010742 353 
010743 042 
130 
022 

010746 041 
060 
024 

010751 011 
010752 176 
010753 022 
010754 043 

r:i|i0tSf^.^023/ 

010760 053 
010761 353 
010762 052 
130 
022 

010765 175 
010766 022 
010767 023 
010770- 174 
010771 022 
010772 052 
134 
022 

010775 353 
010776 041 
060 
034 

01 1001 031 
01 1002 136 
01 1003 043 
01 1004 126 
•01 1005 053 
01 1006 353 


JZ 


NOP 
DAD H 
SHLD BITRN 


XCHG 


MOV 

L/C 

MOV 


DAD 

H 

MOV 

C*L 

MOV 

BaH 

LXI 

H/XREAL 


DAD 

D 

MOV 

E>M 

INX 

H 

MOV 

DvM 

DCX 

XCHG 

"H 

SHLD 

TREAL 

LXI 

H>XREAL 


DAD 

B 

MOV 

A^M 

STAX 

D 

INX 

H 

INX 

D 

MOV 

A^M 

STAX 

D 

DCX 

XCHG 

H 

LHLD 

TREAL 

MOV 

A#L 


STAX D 
I NX D 


MOV A>H 


STAX 

D 

LHLD 

BITRN 


XCHG 

LXI 

H^XIMAG 

DAD 

D 

MOV 


INX ' 

H 

MOV 

D^M 

DCX 

H 

XCHG 




01 1007 042 


SHLD ' 

TIMAG 

132 




022 




01 1 0 1 2 041 


LXI 

H^XIMAG 

060 




034 




011015 011 


DAD 

B 

011016 176 


MOV 

.A^M 

01 1017 022 


STAX 

D 

01 1020 043 


I NX 

H 

011021 023 


INX 

D 

01 1022 176 


MOV 

A/M 




011023 022 STAX D 

01 1024 053 


DCX 

H 

01 1025 353 


XCHG 


011026 052 


LHLD 

TIMAG 

132 




022 




011031 175 


MOV 

A* L 

011032 022 


STAX 

D 

011033 023 


INX 

D 

011034 174 


MOV 

A, K 

011035 022 


STAX 

D 

011036 303 


. JM P 

. LPCON 

054 




022 




011041 175 

CKECAi 

MOV 

A^L 

011042 27'! 


CMP 

C 

01 1043 312 


JZ 

LPCON 

054 




022 




011046 372 


JM 

LPCON 

054 




. 022 




011051 303 


JMP 

EXCHG 

317 




021 




011054 000 

LPGON : 

NOP 


011055 301 


POP 

B 

011056 003 


INX 

B 

011057 052 


LHLD 

N 

012 




016 




011062 174 


MOV 

A#H 

01 1063 270 


CMP 

B 

01 1064 312 


JZ 

SCHEC 

072 




022 




011067 303 


JMP 

XLPl 

301 




021 




01 1072 175 

SCHEC ; 

MOV 

A* L 

01 1073 271 


CMP 

C 

011074 312 


JZ 

XEND 

102 




022 




011077 303 


JMP 

XLPl 

301 




021 




011 102 31 1 

XEND : 

RET 




FftGE 008 


L0G2N CALCULATES POWER OF 2 




NUMBER 

IN 

L) REG. 



POWER 

OF 2 IN 

B REG. 

01 1103 

257 

L0G2N: 

XRA 

A 

01 1 104 

001 


LXI 

B^010B 


010 





000 




01 1 107 

175 


MOV 

A^L 

01 1 1 10 

037 

POWER : 

BAR 


01 1 1 1 1 

332 


JC 

POUT 


127 





022 




011 11 4 

004 


INR 

B 

01 1 1 15 

015 


DCR 

C 

01 1 1 16 

302 


JNZ 

POWER 


110 





022 




01 1 121 

016 


MVI 

C^010B 


010 




01 1 123 

174 


MOV 

A>H 

01 1 124 

303 


JMP 

POWER 


i 10 





022 




01 1 127 

311 

POUT t 

RET 


01 1 130 

000 

TREAL s 

DW 

0B 


000 




01 1 132 

000 

TIMAG 1 

DV 

0B 


000 




011134 

000 

BITRN: 

DW 

0B 


000 






r 


. ' 1 ’ 

- --n 

b 3 


nr 

t-. '' 


- . . 

^ 


«- 


i 

^ . 




y ■ ■• /* 

lu ^ » ♦ 

: 


■ i , V . 


'•"t * ^ - ’T 

1 *■ * 

: 3 

j '2 





^ ' } 

fr/.. ' ' -1: 

' .r ' 


:i :3 

-t 

1 '-"Jht ♦ 

t ;.i 


c ; '-.-S' 

-ni “T| 

'nr 1 


L'.J 

' ■> ■ ■ . 


'*'‘'-0.. "1,' 


. -» 



■ JL* * ... - . * i ♦ 



’ ~3 


INPUT ; 
IF ODD 

ORG 

ROUTINE 

PARITY 

11150B 
READS . 
ENCOUN 

011150 

076 

001 

INPUT: 

MVI 

A^01B 

011152 

323 

010 

- 

OUT 

10B 

01 1 154 

333 

010 

INLP : 

IN 

10B 

011 156 

247 


ANA 

A 

011157 

362 

154 

022 


JP 

INLP 

011162 

333 
01 1 


IN 

1 IB 


37.6 

215 


CPI 

215B 

01 1 166 

312 

150 

022 


J2 

INPUT 

01 1 171 

376 

012 


CPI 

12B 

01 1173 

3*12 

150 

022 


JZ 

INPUT 

01 1 176 

376 

377 


CPI 

377B 

01 1200 

312 

150 

022 


JZ 

INPUT 

01 1203 

376 

254 


CPI 

254B 

011205 

312 

150 

022 


JZ 

INPUT 

01 1210 

376 

000 


CPI 

0B 

011212 

312 

150 


JZ 

INPUT 


JLINE FEED/CONTINUE 




022 





011215 

306 


ADI 

0B 


- 

000 





011217 

342 


JPO 

ERROR 



022 





011222 

346 


ANI 

177B 



177 





011224 

376 


CPI 

7 IB 



071 





011226 

312 


JZ 

NUM ■ 



237 






022 





011231 

332 



NUM 



237 






022 





011234 

326 

, ALPHA J 

SUI 

67B 



067 

? 




01 1236 

31 1 


RET 



01 1237 

326 

NUM : 

SUI 

60B 



060 





01 1241 

311^ 


RET 



01 1242 

166 

ERROR : . 

HLT 


DATA TO ZERO 


initialazation of lsb part of 

01 1243 

001 

ZEROM ; 

LXI 

B>00B 



000 






000 





011246 

076 

ZERO'! * 

MV I 

A> 0 0B 



000 





011250 

167 


MOV 

M^A 


01 1251 

043 

, . ' 

INX 

H 


011252 

043 


INX 

H 



003 


INX 

B 



02 1 


LXI 

D^N 



012 






016 





01 1257 

032 


LDAX 

D 


01 1260 

271 . 


CMP' 

C 


01 1261 

302 


JNZ 

ZEROl 



246 






022 





01 1264 

023 


INX 

D 


01 1265 

032 


LDAX 

D 


01 1266 

270 


CMP 

B 


01 1267 

302 


JNZ 

ZEROl 



246 






022 





01 1272 

311 


RET 


TAPE 

ROUTINE READS DATA FROM PAPER 

01 1273 

043 

READ : 

INX 

H 


01 1274 

000 

READ 1 ; 

NOP 



01 127 5 

001 


LXI 

B,00B 



000 






000 





011300 

315 

READ2 : 

CALL 

INPUT 



150 



022 


DATA IS IN HEX 


011303 

007 


RLC 


011304 

007 


RLC 


011305 

007 


RLC 


011306 

007 


RLC 


011307 

137 


MOV 

'e^a 

01 1310 

315 


CALL 

INPUT 


150 





022 




011313 

263 

* ■ ’ ■ 

ORA 

E 

01 1314 

362 


■' JP 

PLUS 


324 




011317 

022 

137 

MINUS : 

MOV 

E^A ^NEGATIVE NUMBERIN SIGN MAGNITUDE 

01 13'20 

257 


XRA 

^ JCONVERTED TO TWOS COMPLIMEN 

011321 

223 


SUB 

E 

011322 

366 


ORI 

200B 


200 




01 1324 

167 

PLUS : 

MOV 

M^A 

01 1325 

043 


,INX 

H 

01 1326 

043 


I NX 

H 

01 1327 

003 


INX 

B 

011330 

021 


LX I 

D>N 


012 





016 

, , , t 



01 1333 

032 


LDAX 

D 

011334 

27 1 


CMP 

C 

011335 

302 


JNZ 

READ2 


300 

#3( 0 O 





023 


INX 

D 

011341 

032 


LDAX 

D 

01 1342 

270 


CMP 

B 

011343 

302 


JNZ 

READ2 


300 





022 




011346 

311 


RET 


01 1350 



ORG 

1 1350B 



ROUTINES INTFP AND PENT CONVERTS INTEGER 



TO ITS 

REAL 

EQUIVALENT AND PRINT ON TELETYPE 

01 1350 

000 

I NT : 

DW 

0B 


000 




01 1352 

000 

CSTN6 : 

DB 

0B 

01 1353 

076 

INTFP: 

MVI 

A*06B IFORMATTING OF OUTPUT 


006 




01 1355 

062 


STA 

CSTN6 


352 





022 




01 1 360 

046 


MVI 

Hj. 0B 


000 




011362 

072 


LDA 

LL 


016 





016 




011365 

157 


MOV 

L^A 


DAD 

LXI 


H 

D# FLOAT 


011366 051 
011367 021 
277 
023 

011372 031 
011373 136 
011374 043 
011375 126 
011376 353 
011377 042 ■ 

350 

022 

011402 315 
370 

023 

011405 0^1 
060 

024 

011410 315 

025 
023 

011413 315 
370 
023 

011416 041 
060 
034 

011421 315 
§25 

023 . 

01 1424 311 

|s||;|i«05'jt0i ,'PBNT 
000 
000 

011430 136 PLPl 
011431 043 
011432 126 
011433 345 
01 1434 305 
01 1435 052 
350 
022 

011440 315 
121 
016 

01 1443 072 
055 
016 

011446 346 
001 

011450 312 
066 
023 

011453 315 


DAD D 
MOV EjM 
INX H 
MOV D>M 
XCHG 

SHLD I NT 


CALL CBLF ICARRIAGEBETURN* LINE FEED 


LXI H^XREAL 


CALL PRNT 


CALL CRLF 


LXI H^XIMAG 


CALL PRNT 


RET 

LXI B*0B 


MOV E^M 
INX H 
MOV D^M 
PUSH H 
PUSH B 
LHLD I NT 


CALL DPMUL 


LDA FLAG 


ANI 01B 
JZ PLP2 


CALL COMPL 


01 1456 

377 

016 

016 


MVI 

01 1460 

055 

315 


CALL 

01 1463 

355 

023 

303 

» 

JMP 

011466 

073 

023 

016 

PLP2 I 

MVI 

011470 

040 

315 

it 


CALL 

01 1473 

;W ’3' 

023 

257 

PLP2 1 ; 

XRA 

01 1474 

051 


DAD 

01 1475 

051 


DAD 

01 1476 

016 


MVI 

01 1500 

021 

315 


CALL 

011503 

170 

023 

110 


MOV 

01 1504 

315 


CALL 


325 



011507 

023 

112 


MOV 

01 1510 

315 


CALL 


325 

023 

016 


MVI 

01 1515 

056 

315 


CALL 

01 1520 

355 

023 

1 13 


MOV 

01 1521 

315 


CALL 

01 1524 

325 

023 

016 


MVI 

011526 

040 

315 


CALL 

01 1531 

355 

023 

041 


LXI ■ 

011534 

352 

022 

065 


DCR 

011535 

302 


JNZ 

01 1540 

145 

023 

315 


CALL 


370 




C^55B ;for printing minus 
GO 

PLP2 1 

G#40B ; printing BLANK 

A 

H 

H 

C^021B 

BNBCD 

C#B 

OUT 

C^D 

OUT 

C^56B ;PRINT DECIMAL POINT 
CO 

C^E 

OUT 

C>A0B 

CO 

H*CSTN6 

M 

PLP4 

CRLF 




023 



011543 

066 

MVI 

M^06B 


006 



011545 

301 

PLP4 s POP 

B 

01 1546 

341 

POP 

H 

011547 

043 

INX 

H 

011550 

003 

INX 

B 

011551 

021 

LXI 

D>N 


012 




0 16 

* • 


01 1554 

032 

LDAX 

D 

011555 

271 

C,MP 

C 

011556 

302 

JNZ 

PLPl 


030 




023 



01 1561 

023 

INX 

D 

011562 

032 

LDAX 

D 

011563 

270 

CMP 

B 

01 1564 

302 

JNZ 

PLPl 


030 




023 



011567 

311 

RET 



BNBGD CONVERTS 16 BITS BINARY TO BCD 
BINARY NUMBER IN <H,L) REG. 

REG. C SHOULD HAVE IN IT 20 OCTAL 
result IN B>D..E REG. 


01 IS*^® 

345 

BNBCD: 

PUSH 

H 

01 1571 

041 


LXI 

H* TEMPI 


274 





023 




011574 

006 


MVI 

B>03B 


003 




if 1576 

066 

BCDl { 

MVI 

M^00B 


000 




011600 

043 


INX 

H 

011601 

005 


DCR 

B 

011602 



JNZ 

BCDl 


176 





023 




011605 

041 

BCD2 ! 

LXI 

H^ TEMPI 


274 





023 




01 1610 

006 


MVI * 

B>03B 


003 




01 1612 

216 

BCD3 J 

ADC 

M 

01 16 13 

047 


DAA 


01 16 14 

322 


JNC 

BCD4 


224 





023 




01 1617 

026 


MVI 

Dv0lB 


00 1 




01 1621 

303 


JMP 

BCD 5 


226 





023 




01 1624 

026 

BCD4 J- 

MVI 

D>00B 



000 


011626 206 

BCD5 : 

ADD 

M 

011627 047 


DAA 


011630 167 


MOV 

M#A 

011631 322 


JNC 

BCD6 

241 




023 




011634 036 


MVI 

E#01B 

001 




011636 303 


JMP 

BCD7 

243 




023 


: :i'' ■ 


011641 036 

BCD6 : 

MVI 

E^00B 

000 




011643 172 

BCD7 : 

MOV 

A/D 

011644 263 


ORA 

E 

011645 043 


I MX 

H 

011646 005 


DCR 

B 

011647 302 


JNZ 

BCD3 

212 




023 




011652 341 


POP 

H 

011653 257 


XRA 

A 

011654 051 


DAD 

H 

01 1655 345 


PUSH 

H 

011656 015 


DCR 

C 

011657 302 


JNZ 

BCD2 





023 




011662 341 


POP- 

H 


RESIST 

IN TEMPI “TEMP3 



B«D>E 

REG. 

0iT:6Miteii 


, LXI 

H/TEMP3 







MOV 

B/M 



DCX 

H 



MOV 

D/M 

011671 053 


DCX 

H 

011672 136 


MOV 

E/M 

011673 311 


RET 


01 1674 000 

TEMPI : 

DB 

0B 

011675 000 

TEMP2 : 

DB 

0B 

011676 000 

TEMP3 : 

DB 

0B 


FLOAT STORES CONSTANT FOR CONVERTING 

FROM INTEGER TO REAL 

FOR SAMPLE POINTS 0 TO 1024 


01 1677 

000 FLOAT: 
000 

DW 

00B 

01 1701 

062 

DW 

62B 


000 



01 1703 

144 

DW 

144B 


000 



011705 

310 

DW 

310B 


000 


- 

01 1707 

220 

DW 

620B 



011711 Sli 


DW 

1440B 

003 




011713 100 


DW 

3100B 

006 




011715 200 


DW 

6200B 

014 




011717 000 


DW 

I4400B 

031 




011721 000 


DW 

31000B 

062 




011723 000 


DW 

62000B 

144 

OUT 

OUTPUTS 

TWO CHARACTERS 

011725 171 

OUT 

t MOV 

A*C 

011726. 101 


MOV 

B>C 

011727 346 


ANI 

360B 

360 




011731 017 


RRC 


011732 -017 


RRC 


011733 017 


RRC 


011734 017 


RRC 


011735 306 


ADI 

. 60B 

060 




011737 117 


MOV 

CjA 

011740 315 


CALL 

CO 

355 




023 




011743 170 

OUTl 

: m6v 

A^B 

011744 346 


ANI 

17B 

017 




011746 306 


ADI 

60B 

060 




011750 117 


MOV 

C^A 



CALL 

CO 











RET 



CO PRINTS A CHARACTER ON TELE TYPE 
WHEN CALLED CHARACTER BE IN C REG. 

011755 333 CO : IN 12B 

012 

011757 346 ■ ANI 200B 

200 

01 176 1 312 JZ CO 

355 
023 

01 1764 171 MOV A^C 

01 1765 323 OUT 13B 

013 

011767 311 RET 

CRLF FOR LINE SPACING 

011770 016 CRLF : MVI C^15B 

015 

011772 315 
355 


CALL CO 



0117 11 ill DW 

003 

011713 100 D¥ 

006 

011715 200 D¥ 

014 

011717 000 DW 

031 

011721 000 DW 

062 

011723 000 DW 

144 

OUT OUTPUTS 
011725 171 OUT t MOV 

011726, 101 MOV 

011727 346 ANI 

360 

011731 017 RRC 

011732 017 RRC 

011733 017 RRC 

011734 017 RRC 

011735 306 ADI 

060 

011737 117 MOV 

011740 315 CALL 

■ 355 , ■ 

023 

011743 170 OUTl ; MOV 

011744 346 ANI 

017 

011746 306 ADI 

060 

011750 U7 ' MOV 

CALL 

jillifeV'lt'i: ' RET 


CO PRINTS A 
WHEN CALLED 


011755 333 CO : IN 

012 

011757 346 ^ ANI 

200 

01 1761 312 JZ 

355 
023 

011764 171 MOV 

01 1765 323 OUT 

013 

011767 311 RET 

CRLF FOR LINE 

01 1770 016 CRLF ; MVI 

015 

011772 315 GALL 

355 


1440B 

3100B 

6200B 

14400B 

31000B 

62000B 

TWO CHARACTERS 
A,C 
B^C 
360B 


60B 

C^ A 
CO 


A^B 

17B 

60B 

C^A 

CO 


CHARACTER ON TELE TYPE 
CHARACTER BE IN C REG. 
12B 

200B 

CO 


A#G 

13B 


SPACING 

C^15B 

CO 



MV I 


Cj 12B 


023 

011775 016 
012 

011777 315 CALL CO 

355 

023 

012002 311 RET 

012010 ORG 12010B 

MAIN PROGRAM STARTS 

012010061 LXI SP>27777B ;INIALIZEZTACK POINTER 

377 
057 

? number OF sample POINTS ENTERED IN (H^L)REG 
012013 041 H#200B 

200 
000 

FOR 128 SAMPLE POINTS 
LOWER BYTE IN LOCATION 12014 
HIGHER BYTE IN LOCATION 12015 
012016 042 SHLD N 

012 

016 ■ 

012021 041 LXI H/XREAL 

060 

024 

012024 315 CALL ZEROM IFOR MAKING LSB PART OF DATA ZERO 

■ 243 

' 022 ■ ■ ' ■ 

012027 041 LXI H/XIMAG 

' 060 - " ' : 

, ;.:' 034 / 

012032315 CALL ZEROM 

243 

060 
. 024 

012040 315 
273 

622 

012043 041 LXI H>XIMAG 

060 
034 

012046 315 CALL READ 

|^®E 2 X 02 S 

273 

022 

012051 315 CALL FFTIN IFOR FFT 

000 
014 

012054 315 CALL INTFP IFOR INTEGER TO REAL CONVERSION 


LXI H^XREAL 

CALL READ IFOR READING REAL DATA FROM TAPE 



353 

022 

AND PRINTING ON TELE TYPE 
012057 166 HLT 

RESERVE LOCATIONS FOR 1024 SAMPLE POINTS 
012060 XREAL: DS 4000B 

016060 XIMAG: DS 4000B 

END 

NO OF ERRORS! 00 



A 55283 



This book Is to be returned on the 
date last stamped* 








